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ABSTRACT 


A  combined  experimental/analytical  research  program  aimed  at 
extending  the  static  J-lntegral  to  the  case  of  dynaml''  Inelastic 
fracture  Is  described.  This  program  was  a  first  step  .  '.ward  the 
goal  of  quantifying  the  fracture  response  of  naval  structures 
containing  flaws  subjected  to  dynamic  loading.  One-inch  thick 
HY-80  baseplate  was  used  because  of  Its  known  toughness  and  its 
widespread  application  In  naval  structures. 

For  the  experimental  work,  a  drop  weight  test  fixture  was  de¬ 
signed  that  permitted  controllable,  high  rate  loading  of  precracked 
compact  tension  specimens.  A  series  of  specimens  was  tested,  and 
applied  load,  load  point  displacement,  and  back  face  strains  were 
routinely  measured  with  improving  confidence  and  accuracy  as  the 
program  progressed.  The  final  measure  required  to  obtain  a  critical 
value  of  the  J-integral  was  the  time,  load  or  load  point  displacement 
at  which  fracture  Initiated.  This  proved  to  be  the  most  difficult 
event  to  measure.  After  much  trial  and  error,  a  multiple  specimen 
approach  was  adopted  in  which  specimens  were  loaded  up  to  and  just 
beyond  the  point  of  crack  initiation.  Using  this  method,  an  apparent 
dynamic  Jjc  HY-80  steel  was  "bracketed”.  The  resulting  value 
was  approximately  1000  In.-lb/in.^  versus  an  average  static 
of  1110  in.-lb/in.2  for  the  same  plate. 

For  the  analytical  work,  two-dimensional  finite  element  analyses 
were  performed  in  conjunction  with  the  experiments  In  an  effort  both 
to  further  develop  analytical  techniques  for  dynamic  fracture  problems 
and  to  provide  information  on  the  specimen  response  in  this  test 
program.  The  complex  state  of  stress  and  the  presence  of  side  grooves 
in  the  specimens  complicated  the  application  of  two-dimensional 
analyses.  These  problems  were  overcome,  but  the  precision  of  the 
analyses  performed  remains  uncertain  pending  further  verification, 
perhaps  using  three-dimensional  analysis  methods. 

The  analyses  performed  displayed  variable  correlation  with  the 
experimental  results.  This  variability  was  dependent  on  the  yield 
stress  input  to  the  analysis.  Surprisingly  large  differences  in 
load  point  displacement  resulted  from  small  changes  in  yield  stress. 

The  analyses  also  showed  the  load  point  displacement  to  be  very  sensi¬ 
tive  to  the  modeled  crack  depth.  Changes  of  less  than  1%  in  the 
ratio  of  crack  length  to  specimen  width  resulted  in  load  point  dis¬ 
placement  variations  of  roughly  40%.  This  sensitivity  to  stress-strain 
properties  and  crack  depth  was  thought  to  be  largely  due  to  the  deeply 
cracked  nature  of  the  specimens  in  this  program,  but  it  remains  a 
point  of  concern  and  interest  due  to  the  potential  effect  on  similar 
fracture  analyses. 

ADMINISTRATIVE  INFORMATION 

This  project  was  initiated  and  supported  primarily  by  the  David  W.  Taylor  Naval 
Ship  Research  and  Development  Center  (DTNSRDC)  Independent  Research  Program,  Element 
hllS2N,  Task  Area  ZR  023  0301,  Work  Unit  1720-410. 
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INTRODUCTION 


The  initiation  and  growth  of  cracks  in  critical  naval  structures  has  always 
been  a  concern  to  those  who  design  scantlings  and  select  materials  and  fabrication 
piocedures.  Many  designs,  however,  notably  lack  quantitative  consideration  of  de¬ 
fects  ("How  harmful  Is  a  given  flaw  under  given  loading?").  Instead,  reliance  has 
been  placed  on  very  tough  (defect  tolerant)  materials,  high  quality  welding,  rigid 
weld  inspection,  and  follow-on  weld  surveillance.  This  state  of  affairs,  and 
ongoing  research  to  improve  the  situation,  have  been  discussed  by  a  number  of 

Navy  researchers . 

Cracklike  defects  in  tough  materials  have  not  been  addressed  quantitatively 
because  adequate  technology  simply  does  not  yet  exist.  Beca.se  the  materials  in 
question  are  extremely  tough,  the  relatively  well-established  principles  of  linear 
elastic  fracture  mechanics  are  not  often  applicable  except  In  fatigue  analysis. 

The  emerging  science  of  nonlinear  fracture  mechanics,  on  the  other  hand,  shows 
great  promise  for  dealing  quantitatively  with  cracks  in  tough  materials,  that  is, 
for  assessing  the  margin  of  safety  associated  with  a  given  flaw  for  a  given  static 
or  dynamic  load. 

Within  the  past  decade,  research  In  nonlinear  fracture  mechanics  has  con¬ 
centrated  largely  on  static  loading,  with  the  so-called  "J-i ntegral becoming 
widely  accepted  as  a  fracture  toughness  measure.  Standard  procedures^  were 
developed  for  J-lntegral  testing  under  static  conditions,  and  such  tests  have  been 
conducted  routinely  at  this  Center  for  some  time.®  The  early  rationale  behind 
this  work  was  the  development  of  methods  for  fracture  toughness  testing  of  more 
ductile  metals,  which  are  beyond  the  linear  elastic  (plane  strain)  constraints 


*References  are  listed  on  page  121  . 


imposed  by  ASTM  Standard  E399.^  Subsequently,  it  was  recognized^O  that  the 
j- 1 ntegral  also  could  serve  as  a  fracture  initiation  criterion  in  real  cracked 
structures  (as  opposed  to  standard  test  specimens),  provided  one  could  calculate 
its  value  as  a  function  of  applied  load. 

The  loading  of  most  interest  for  many  naval  structures  is  dynamic,  however, 
and  this  area  of  nonlinear  fracture  mechanics  has  been  subject  of  practically  no 
research  other  than  the  recent  work  described  in  References  11  through  14.  This 
report  describes  a  basic  experimental  and  analytical  research  program  aimed  at 
extending  the  J-integral  approach  to  the  analysis  of  inelastic  fracture  initia¬ 
tion  under  dynamic  loading.  Our  research  was  motivated  by  the  following 
questions : 

1 .  Can  the  J-integral  concept  be  used  to  characterize  inelastic 
fracture  initiation  under  dynamic  loading  conditions? 

2.  How  can  the  critical  value  of  J  (Jic)  be  measured  under  dynamic 
loading  conditions? 

3.  Is  Jjc  dependent  on  loading  rate  for  steels  typically  used  in  sub¬ 
marine  hull  construction? 

Such  questions,  it  was  found,  are  not  easy  to  answer.  This  report  covers  "trial 
and  error"  research  conducted  over  4  years.  Considerable  detail  of  the  experi¬ 
mental  and  analytical  programs  is  included  so  that  future  researchers  may 
benefit  from  our  experiences,  both  good  and  bad. 

BACKGROUND 

The  J-integral  is  used  in  nonlinear  fracture  mechanics  to  measure  the 
intensity  of  the  stress  and  strain  fields  near  the  tip  of  a  crack  that  tend  to 
drive  the  crack  forward.  As  loading  on  a  cracked  structure  is  increased,  the 


J-integral  also  increases  (usually  In  a  nonlinear  fashion).  When  J  reaches  Its 
critical  value,  Jic»  crack  growth  will  initiate.  With  a  further  increase  in 
load,  stable  crack  growth  or  unstable  crack  propagation  (complete  fracture) 
will  occur,  depending  on  the  method  of  loading  (displacement  control  versus 
load  control)  and  on  the  structural  compliance. 

The  J-integral  is  thus  not  a  failure  criterion.  If  the  structure  and 
loading  cause  crack  growth  to  decrease  the  crack  driving  force  (a  common  situ¬ 
ation),  then  failure  will  not  occur  when  J  reaches  Jjc* 

Physically,  the  J-lntegral  is  defined  as  the  rate  at  which  total  potential 
energy  of  a  cracked  body  changes  with  respect  to  developing  crack  surface  area, 
and  is  a  nonlinear  generalization  of  the  elastic  strain  energy  release  rate  G.  It 
is  defined  for  two-dimensional  through-cracked  structures  as 

j  =  -ilU  (1) 

B  9a 

where  B  is  the  thickness  of  the  structure,  U  is  the  potential  energy,  and  a  is  the 
crack  depth.  As  shown  in  Figure  1  for  compact  tension  specimens  of  unit  thickness, 
jAa  at  the  load  level  Pj  is  simply  the  shaded  area  between  the  load-displacement 
curves  for  two  specimens  having  crack  lengths  of  a  and  a  +  Aa.l5 

In  addition  to  this  experimental  interpretation,  the  J-integral  has  a  mathe¬ 
matical  definition  as  given  by  Rlce^  for  two-dimensional  bodies  containing  a  mode 
I  through  crack; 


J  =  y  Wdy  -  ds  (i  =  x,y) 


Here  F  is  any  closed  path  surrounding  the  crack  tip,  W  is  the  strain  energy  density 


(■■f" 


Ij  are  the  tractions  acting  on  the  contour  (Tj  =  n  j ) , 
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COMPACT  TENSION  SPECIMEN  (UNIT  THICKNESS) 


Figure  1  -  Physical  Interpretation  of  the 
J-Integral  for  a  Compact 
Tension  Specimen 

p-ij  are  strains,  are  stresses,  nj  are  the  components  of  the  unit  normal  to  the 
path,  u^^  are  displacements,  and  x,  y,  and  ds  are  as  defined  in  Figure  2.  Rice  showed 
that  for  linear  elastic  bodies,  and  for  bodies  to  which  the  deformation  theory  of 
plasticity  applies,  the  J-lntegral  is  independent  of  the  path  chosen  for  its  evalua¬ 
tion.  For  the  more  realistic  Incremental  theories  of  plasticity,  the  J-integral 
is  practically  path-independent  so  long  as  the  loading  increases  monotonically  and 
remains  nearly  proportional.  Thus  the  J-integral  is  an  ideal  parameter  for  calcula¬ 
ting  crack  driving  force,  since  the  researcher  may  choose  a  path  that  is  remote  from 
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Figure  2  -  The  Arbitrary  Path  F  About  a  Crack  Tip 


the  crack  tip,  where  accurate  stress  analysis  is  easier,  rather  than  analyzing  the 
vicinity  of  the  crack  tip,  which  is  difficult.  In  general,  nonlinear  finite  element 
analysis  is  used  to  calculate  the  J-integral  as  a  function  of  increasing  load  in 
realistic  structures  (as  opposed  to  test  specimens).  As  an  example,  the  DTNSRDC 
PAPST  (plastic  Axisymmetric/^lanar  STructures)  finite  element  computer  program 
has  a  built-in  capability  to  calculate  the  J-integral  on  up  to  ten  user-specified 
paths  about  a  crack  tip. 

For  compact  tension  specimens  like  that  shown  in  Figure  1 ,  the  experimental 
value  of  the  J-integral  may  be  related  to  the  load  versus  load  point  displacement 
curve  in  the  manner  shown  by  Merkle  and  Corten^^; 


J  = 


2(l+a) 


B(W-a)(l+a 


^  ,  2a(l-2.a-a^  ^ 

■)  J  y  B(W-a)(l+a^)^  i  ^ 


(3) 


.V. 


-1, 


V. 


In  the  above,  /pduy  is  the  area  under  the  load-displacement  curve  (work), 
/ UydP  is  the  complementary  work,  and  a  is  given  by 

v2  ,,  -11/2 


a  = 


(2a_\ 

Vw-a  j  V-a. 


2a 


W  -  a 


-  1 


(4) 


where  W  and  a  refer  to  the  specimen  dimensions  shown  in  Figure  1. 

Under  slow  test,  the  critical  value  of  J  may  be  determined  in  a  single  specimen 
test  by  noting  the  load  or  load  point  displacement  at  which  fracture  begins.  This 
is  usually  done  by  the  unloading  compliance  method  described  by  Joyce  and  Gudas,® 
wherein  the  load  is  decreased  by  about  lOX  periodically  during  a  test  to  check  for 
a  change  in  elastic  specimen  compliance  that  indicates  fracture  initiation.  The 
load-displacement  record,  together  with  Equation  (3),  then  gives  Jxc*  fact, 
for  deeply  cracked  specimens.  Equation  (3)  can  be  reduced  tol® 


1  +  a 


2A„ 


1  +  a 


2  B(W-a) 


(5) 


where  Ax  is  the  total  area  under  the  load-load  point  displacement  curve  diagrammed 
in  Figure  3. 

To  summarize,  we  have  outlined  two  uses  of  the  J-integral  in  the  static  realm: 
(1)  as  an  experimental  parameter  (beyond  linear  elastic  fracture  mechanics)  for 
measuring  the  toughness  of  materials  loaded  into  the  elastic-plastic  range,  and  (2) 
as  a  mathematical  parameter  that  can  be  calculated  for  realistic  elastic-plastic 
cracked  structures.  This  duality  affords  a  capability  for  design  against  fracture 
initiation  by  comparing  calculated  values  of  J  against  the  experimentally  measured 
critical  value  of  the  J-integral,  Jjc*  Until  now,  however,  this  idea  has  been 
Largely  confined  to  the  static  case. 
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Figure  3  -  Calculation  of  J  from  a  Load- 
Displacement  Record 


In  the  work  described  in  this  report  we  sought  to  extend  the  J-lntegral  concept 
to  dynamic  inelastic  fracture  of  large-scale  naval  structures.  The  specific  objec¬ 
tives  of  this  research  were  as  follows: 

1 .  Design  an  experimental  technique  by  which  standard  compact  tension  Jjq 
specimens  may  be  rapidly  loaded  at  various  rates. 

2.  Develop  Instrumentation  techniques  capable  of  measuring  fracture  phenomena 
that  occur  in  less  than  a  millisecond.  Specifically,  load  versus  time  and  load 
point  displacement  versus  time  must  be  accurately  measured  in  order  to  use  equations 
like  (4)  and  (5)  to  construct  the  curve  of  the  J-integral  versus  time. 

3.  Develop  a  technique  for  determining  when  fracture  begins  in  terms  of  time, 
load,  or  displacement  in  order  to  estimate  the  value  of  J-[c  under  dynamic 
conditions , 
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4.  Validate  analytical  techniques  (nonlinear  dynamic  finite  element  fracture 
analyses)  to  assist  in  guiding  the  experimental  program  and  to  make  possible  the 
calculation  of  J  In  realistic  cracked  ship  structures. 

5.  Determine  the  difference,  if  any,  between  the  values  of  Jjq  under  static 
and  dynamic  loading  for  typical  (HY)  steel  structural  materials. 

Our  approach  and  progress  in  achieving  these  goals  are  described  in  this  report. 

EXPERIMENTAL  TEST  PROGRAM 

TEST  SPECIMEN 

The  J-integral  has  been  Identified  as  an  appropriate  criterion  for  fracture 
initiation  for  many  tough  materials  under  static  and  quasi-static  load  conditions.^® 
Oar  task  was  to  determine  this  same  property  for  dynamic  load  conditions  at  rates 
typical  of  those  associated  with  an  impact  or  an  explosion  attack  on  naval  struc¬ 
tures.  The  deep  notch  compact  tension  (CT)  specimen  used  in  these  studies  was 
similar  to  the  typ:  widely  used  to  study  material  toughness,  perhaps  most  notably 
the  plane-strain  fracture  toughness  properties  as  outlined  in  Reference  9.  Basic 
dimensions  for  the  specimen  used  in  this  work  are  shown  in  Figure  4,  and  manufac¬ 
turing  drawings  are  identified  in  Reference  19. 

Alterations  were  made  to  the  specimen  during  the  test  program.  The  alterations 
Included  modifying  the  holes  to  accept  the  Instrumentation  and  machining  side 
grooves  to  reduce  undesirable  shear  lips  that  tend  to  develop  on  the  specimen  edges. 
A  fatigue  precrack  had  to  be  introduced  into  the  specimen  before  the  test.  The 
controlled  procedure  outlined  in  Reference  9  was  used  to  assure  that  a  crack  of 
proper  depth  and  with  a  sharp  tip  was  obtained  so  that  the  specimen  would  perform 
consistently  and  give  valid  data. 


TEST  APPARATUS 


After  considering  the  various  test  conditions  desired  for  this  study,  we  decid¬ 
ed  that  a  drop  test  technique  would  be  most  appropriate.  The  specimen  would  be  in¬ 
stalled  in  a  stationary  fixture  which  a  free-falling  weight  would  strike  thus  load¬ 
ing  the  specimen  in  the  opening  mode.  This  approach  was  simple  and  offered  many 
Loading  rate  options.  The  applied  load  could  be  varied  by  choice  of  drop  height  and 
the  size  of  the  drop  weight.  Drop  weights  of  115,  60,  and  30  lb  were  available. 

The  heaviest  weight  was  made  of  steel  and  the  two  lighter  ones  of  aluminum.  They 
were  cylindrical,  8  in.  in  diameter,  and  were  constructed  with  a  longitudinal  groove 
on  each  side  to  mate  with  a  vertical  guide  track.  An  impact  buffer  could  be  placed 
between  the  fixture  and  the  impacting  weight.  These  buffers  could  be  designed  with 
different  properties  to  act  as  load  programmers  providing  even  wider  test 
variations . 

Some  preliminary  work  was  done  to  determine  the  extent  to  which  the  loading 
could  be  controlled  by  varying  certain  test  parameters,  namely  the  impact  device  and 
the  drop  height.  Several  different  materials,  including  nylon,  adiprene  (a  rubber 
type  material),  and  lead,  were  investigated.  These  materials  were  configured  into 
various  cross-sectional  areas  and  thicknesses  to  give  different  stiffness  character¬ 
istics.  To  investigate  their  dynamic  characteristics  these  devices  were  then  placed 
on  a  rigid  steel  base  and  were  struck  by  a  60-lb  weight  dropped  from  two  discrete 
heights.  This  facility  was  designed  so  that  it  could  later  be  used  in  the  test 
program.  To  determine  the  performances  of  these  different  arrangements,  the  decel¬ 
eration  pulse  was  measured  at  Impact  with  an  accelerometer  positioned  on  top  of  the 
drop  weight.  These  pulses  are  qualitative  indications  of  the  shape  and  duration  of 
the  generated  force.  They  do  not  truly  represent  the  actual  test  pulses,  because 


the  dynamic  characteristics  of  both  the  test  fixture  and  of  the  specimen  would  alter 


the  load  pulses  obtained.  The  shapes  and  the  range  of  magnitudes  and  durations  for 
the  load  pulses  are  shown  in  Figure  5. 

An  extremely  rigid  test  fixture,  which  permitted  the  specimen  to  be  loaded  in 
an  opening  mode  at  Impact,  was  designed  and  built.  It  was  designed  with  tight 
tolerances  so  that  essentially  only  an  axial  load  would  be  applied.  Special  clevis 
pins  to  engage  the  specimen  were  used  to  eliminate  clearance  and  prevent  extraneous 
signals  from  occurring  on  the  test  measurements.  Figure  6  is  a  sketch  of  this  drop 
test  fixture;  Reference  20  identifies  drawings  used  to  fabricate  the  device. 

INSTRUMENTATION 

The  J-integral  is  a  measure  of  the  intensity  of  stress  and  strain  fields  near  a 
crack  tip  in  a  material  at  a  given  value  of  applied  load.  At  some  load,  a  critical 
condition  known  as  Ojc  is  reached,  where  crack  growth  initiates.  The  J-integral 
can  be  conveniently  obtained  from  a  curve  relating  the  applied  load  and  the  corre¬ 
sponding  load  point  displacement.  To  determine  Jic*  crack  initiation  must  be 
determined  as  a  function  of  time,  displacement,  or  load. 

Special  instrumentation  and  techniques  had  to  be  developed  to  obtain  these 
measurements  because  the  tests  were  to  be  conducted  under  rapid  load  conditions.  A 
large  effort  was  Involved  with  this  phase,  and  frequently,  certain  procedures  were 
found  unsuccessful  and  alternate  approaches  were  needed.  The  load  point  displace¬ 
ment  and  the  crack  Initiation  were  particularly  difficult  to  measure.  A  technique 
was  finally  developed  whereby  this  displacement  could  be  confidently  measured,  but 
the  problem  of  detecting  crack  initiation  in  a  single-specimen  test  was  not  solved 
satisfactorily.  Without  this  capability,  the  alternative  was  to  run  a  series  of 
tests,  with  increasingly  severe  loading,  and  to  determine  the  threshold  of  fracture. 
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It'e  corresponding  J-integral  was  Chen  assumed  to  be  the  critical  value,  Jxc*  This 
approach  was  used  for  two  basically  different  load  histories.  At  this  writing,  26 
HY-80  steel  specimens  have  been  tested  and  42  separate  drops  made.  Table  1  lists 
these  specimens  along  with  the  essential  parameters  (drop  height,  drop  weight,  and 
impact  device)  for  each  test. 

Following  are  discussions  of  the  procedures  and  techniques  used  to  obtain  the 
essential  measurements.  Some  procedures  were  quite  successful  from  the  start  and 
remained  essentially  unchanged  throughout  the  test  series.  Obtaining  some  of  the 
more  difficult  measurements,  however,  required  an  evolutionary  process  of  varying 
instrumentation  and  associated  techniques.  In  the  Interest  of  clarity,  the  next 
sections  describe  only  the  most  accurate  and  reliable  of  the  measurement  methods 
attempted.  Detailed  discussion  of  all  the  techniques  and  instrumentation  used  and 
the  trials  associated  with  obtaining  valid  data  from  them  are  provided  In  the 
appendix. 

Applied  Load 

Figure  7  is  a  sketch  of  the  test  fixture  and  compact  tension  specimen.  As 
shown,  the  load  cell  was  located  on  the  lower  clevis  stem  and  measured  the  force 
transmitted  to  the  specimen.  This  cell  consists  of  a  strain  gage  bonded  to  each 
side  of  the  stem.  They  operate  as  two  active  arms  of  a  four  arm  Wheatstone  bridge 
and  are  connected  in  opposite  arras  of  the  bridge  to  effectively  cancel  any  bending 
components  on  the  signal.  The  cross  section  of  the  clevis  where  the  gages  v/ere 
attached  had  been  reduced  in  order  to  increase  the  sensitivity  of  the  bridge.  Thi 
load  cell  was  calibrated  statically  for  loads  up  to  15,000  lb  and  the  results  are 
shown  in  Figure  8.  Also  shown  in  the  figure  is  the  projected  sensitivity  based  on 
the  cross  section  area  of  the  stem;  note  the  close  correspondence  between  the  two. 


*  See  note  at  end  of  table 


TABLE  1  (Continued) 


Specimen 

No. 

Drop 

No. 

Drop 

We ight 
(lb) 

Drop 

Height 

(in.) 

Element* 

18 

3 

60 

12 

3 

4 

60 

8 

1 

3 

5 

60 

8 

3 

19 

1 

60 

8 

3 

2 

60 

8 

3 

3 

60 

8 

3 

4 

60 

8 

3 

21 

1 

60 

8 

3 

2 

60 

8 

3 

20 

1 

60 

8 

3 

2 

60 

8 

3 

22 

1 

60 

8 

3 

23 

1 

60 

8 

18 

26 

1 

60 

12 

18 

27 

1 

60 

10 

18 

31 

1 

60 

10 

18 

32 

1 

60 

12 

18 

33 

1 

60 

11 

18 

29 

1 

60 

11 

18 

*  Welght/f ixture  Impact  device:  Type  8,  3  in.  dia.  x  in.  thick  Adiprene; 
Type  3,  l^s  in.  dia.  x  l  in.  thick  Nylon;  Type  18,  1  in.  x  1  in.  x  1/8  in. 
thick  lead. 
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Figure  7  -  Initial  Instrumentation  for  Obtaining  Jjc 
(Load  Cell  and  Displacements  D1  and  D2  or  D2') 


This  static  sensitivity  was  used  to  determine  the  force  applied  to  the  specimen 
during  the  dynamic  tests. 

Load  Point  Displacement 

The  second  essential  measurement  was  the  load  point  displacement,  LPD — also 
called  the  load  line  or  the  crack  opening  displacement  (LLD  and  COD,  respectively). 
This  is  measurement  D1  in  the  schematic  of  Figure  7.  Primary  considerations  for 
selecting  this  sensor  were  frequency  response,  displacement  range,  and  size  and  ease 
of  adapting  to  the  test  setup.  A  clip  gage,  which  had  been  used  in  the  past  for 
this  application  under  static  loading,  did  not  offer  these  desired  features. 
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Figure  8  -  Calibration  of  Test  Fixture  Load  C^ll 


After  attempts  with  a  non-contact  eddy-current  displacement  measuring  system 
proved  less  than  successful  due  to  problems  with  the  linearity  of  the  sensor,  a 
system  using  fiber  optic  transmitted  and  reflected  light  to  determine  load  line 
displacement  was  used.  Detailed  descriptions  of  both  of  the  systems  and  the  cali¬ 
bration  and  test  procedures  tried  for  each  are  presented  in  the  appendix. 

The  fiber  optic  system  provided  the  most  consistent  results;  the  displacement 
measurements  obtained  using  it  were  accurate  to  within  10%.  This  system,  known  as  a 
Fotonlc”  sensor,  was  manufactured  by  Mechanical  Technology,  Inc.  The  sensing  probe 


contains  two  sets  of  fiber  optic  filaments  arranged  in  concentric  circles.  One  set 
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of  filaments  transmits  light  from  the  instrument  to  a  target  on  the  lower  face  of 
the  specimen  notch.  The  second  set  receives  the  reflected  light  and  transmits  it 
back  to  the  sensor  unit.  There,  the  reflected  light  is  converted  to  an  electical 
signal  that  is  proportional  to  the  gap  between  the  probe  tip  and  reflecting  surface. 

A  large  effort  was  required  to  obtain  a  good  calibration  curve  for  the  Fotonic 
sensor.  The  primary  problem  was  that,  because  of  sensor  characteristics,  the  effec¬ 
tive  hinge  length^l  of  the  compact  tension  specimen  had  to  be  considered  during  the 
micrometer  calibration  of  the  sensor.  The  opening  of  the  notch  faces  produced  a 
nonparallel  reflecting  surface  which  reduced  the  accuracy  of  the  sensor.  This  was 
accounted  for  by  using  a  steel  rotating  face  during  the  calibration  of  the  sensor  so 
that  the  sensitivity  curve  obtained  was  representative  of  the  test  conditions.  A 
complete  description  of  this  procedure  is  provided  in  the  appendix. 

Detection  of  Crack  Initiation 

The  final  measurement  required  to  determine  Jjc  is  an  indication  in  real  time 
of  crack  initiation  in  the  specimen.  Several  techniques  were  tried  but  none  was 
successful  enough  to  give  reliable  data.  As  a  consequence,  part  of  the  test  program 
emphasized  loading  the  specimen  just  to  the  threshold  of  crack  initiation.  Thus, 
the  measured  J-integral  would  correspond  roughly  to  the  desired  Jic* 

Though  not  finally  successful,  an  extensive  effort  was  made  to  determine  pre¬ 
cisely  when  the  crack  extension  began.  Many  of  the  techniques  tried  and  types  of 
instrumentation  used  showed  promise  for  application  in  similar  work.  The  effort  is 
described  completely  in  the  appendix.  The  various  techniques  used  and  problems  en¬ 
countered  in  this  work  are  described  briefly  below. 

Initially,  an  eddy  current  displacement  sensor  like  that  used  for  the  first 
few  load  point  displacement  measurements  was  placed  in  a  hole  in  the  ligament 
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section  of  the  specimen  just  beyond  the  crack  tip.  Subsequently,  the  sensor  was 
moved  to  the  notch  area  of  the  specimen  just  before  the  crack  tip.  These  locations 
are  shown  as  D2  and  D2'  in  Figure  7.  We  had  hoped  that  these  sensors  would  indicate 
crack  initiation  by  showing  a  significant  change  in  the  slope  of  the  ligament  or 
notch  displacement  plotted  versus  the  load  point  displacement  response.  No  trend 
indicating  crack  initiation  was  determined  because  the  sensor  in  the  ligament  was 
Insensitive  to  crack  extension  until  the  crack  reached  the  hole  containin?"  the 
sensor.  Also,  the  notch-mounted  sensor  displayed  no  marked  change  in  slope  when 
plotted  versus  the  load  point  displacement  sensor. 

The  failure  of  the  displacement  response  method  shifted  the  emphasis  to  the  use 
of  more  direct  methods  at  or  near  the  crack  tip  to  determine  initiation.  Three 
different  procedures  were  employed:  crack  detection  strain  gages,  crack  propagation 
gages  and  an  optical  moire  fringe  method. 

During  the  test  program,  we  observed  substantial  plasticity  at  the  side  faces 
of  the  specimens.  There  was  concern  that  this  could  Influence  the  fracture  charac¬ 
teristics.  To  eliminate  this  effect,  V-grooves  were  machined  on  both  sides  of  the 
specimen.  Finally,  to  accommodate  bondable  gages,  wide  grooves  with  a  radius  at  the 
bottom  were  used,  although  this  arrangement  perhaps  was  not  as  successful  in  elimi¬ 
nating  the  formation  of  shear  lips  upon  fracture  as  the  V-grooves.  The  configura¬ 
tion  of  the  grooves  is  shown  in  Figure  9. 

Figure  10  shows  the  arrangement  of  the  crack  detection  and  crack  propagation 
gages.  The  crack  detection  gages  were  standard  strain  gages,  applied  to  the  speci¬ 
men  sides  about  1/8  in.  away  from  the  visible  crack  tip.  We  hoped  that  the  plastic 
strains  measured  near  the  crack  tip  would  decrease  markedly  when  crack  extension 
occurred  and  that  these  gages  would  show  this  effect.  Unfortunately,  the  gages 
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Figure  10a  -  Crack  Propagation  Gage 
(Type  MM  TK-09-CPB02-005) 


Figure  10b  -  Strain  Gage  on  Flat  Surface 
(Type  BLH  FAE-03-12-56  ET) 


SIDE  GROOVE 


Figure  10c  -  Strain  Gage  In  Wide  Groove 
(Tvoe  BLH  FAE-03-12-56  ET) 


Figure  lOd  -  "Krak"  Gage  In  Wide  Groove 
(ITT  Type  KG- AOS  and  KG-AIO) 


Figure  10  -  Procedures  Using  Bondable  Gages  to  Detect  Crack  Propagation 

tended  to  come  unbonded  due  to  the  high  strains  occurring  where  they  were  located. 
This  unbonding  produced  the  same  type  of  signal  as  expected  for  the  strain  release 
caused  by  the  extension  of  the  crack.  Because  of  this  unbonding,  consistently 
detecting  crack  initiation  was  Impossible. 

The  crack  propagation  gages  also  suffered  from  this  unbondlng  problem.  These 
gages  are  designed  to  fracture  with  the  crack  extension,  thus  producing  a  measurable 
change  In  gage  resistance.  The  unbondlng  prevented  the  gage  from  following  the 


crack  extension  In  the  specimen. 


High  speed  photography  and  moire  fringes  were  used  briefly  with  the  hope  of 
determining  strains  in  the  specimen  and  crack  initiation.  Because  the  event  was  so 
fast,  the  films  obtained  did  not  have  good  resolution.  In  addition,  the  overall 
plasticity  due  to  the  hinging  of  the  specimen  tended  to  mask  the  moire  fringe 
measurements  in  the  crack  tip  region.  No  conclusive  strains  or  indications  of  exten¬ 
sion  could  be  observed  in  the  crack  tip  region  using  this  technique.  The  use  of 
side  grooved  specimens  precluded  any  further  use  of  the  moire  fringes,  but  high 
speed  films  were  still  selectively  taken  and  compared  with  the  load  point  displace¬ 
ment  measurements. 

The  final  method  used  in  the  effort  to  detect  the  onset  of  crack  extension  was 
direct  current  potential  drop  (DCPD  or  PD).  This  technique  relies  upon  the  increase 
in  electrical  resistance  of  the  cross-sectional  area  of  the  specimen  as  the  crack 
grows  into  the  ligament.  Considerable  work  has  been  done  applying  this  method  in  the 
static  testing  of  compact  tension  specimens.  That  work  was  heavily  consulted  about 
the  intricacies  of  the  electrical  probe  locations  on  the  specimen, 22  and  the  diffi¬ 
culties  associated  with  using  DCPD  during  fatigue  precracking  and  on  already  pre¬ 
cracked  specimens .2 ) >24  attachment  points  for  the  input  current  and  potential 

probes  are  shown  in  Figure  11.  To  verify  the  system,  the  potential  for  various 
crack  depths  (determined  by  compliance)  was  checked  during  precracking  against  cali¬ 
bration  work  on  a  saw  cut  CT  specimen  like  that  described  in  Reference  24.  This 
comparison  is  described,  along  with  a  thorough  description  of  the  potential  drop 
method  as  applied  in  this  experimental  program,  in  the  appendix. 

Problems  arose  due  to  the  variable  effect  that  closure  of  the  fatigue  precrack 
had  on  the  specimen  potential  when  the  specimen  was  unloaded.  We  decided  that  the 
problem  could  be  avoided  by  measuring  the  final  potential  under  load  during  precrack- 
Ing  and  determining  crack  initiation  when  that  value  was  exceeded  during  the  test. 


POTENTIAL  LEADS 


Figure  11  -  Potential  and  Current  Lead 
Attachment  Locations 


When  this  procedure  was  applied  during  the  tests,  however,  totally  unexpected  volt¬ 
age  outputs  were  obtained.  The  voltage  changes  were  much  greater  than  any  of  those 
obtained  during  the  specimen  precracking.  Apparently,  the  instrumentation  added  for 
the  test  (the  D2'  sensor  and  the  common  grounding  of  it)  caused  these  surprising 
results . 

Only  two  specimens  were  tested  using  the  DCPD  procedure  and  both  produced 
extreme  2-  to  3-V  negative  to  positive  potential  changes.  Although  these  results 
were  certainly  incorrect  and  possibly  caused  by  an  electrical  problem  in  the  instru¬ 
mentation,  the  DCPD  method  showed  a  good  deal  of  promise.  Time  and  cost  constraints 
prevented  further  examination,  however. 


Thus,  the  trial  and  error  testing  produced  only  one  measure  of  crack  lnltiatl'*n 
that  would  be  considered  valid:  loading  to  the  threshold  of  crack  initiation  as 
determined  by  breaking  open  the  specimen  after  the  test. 

Summary  of  Instrumentation  Design 

A  viable  instrumentation  plan  was  developed  where  dynamic  J-lntegral  tests 
could  be  conducted  for  different  load  rates.  Techniques  were  developed  to  obtain 
all  essential  measurements  for  the  J-integral  except  crack  initiation.  In  this 
effort,  various  techniques  were  tried,  including  a  secondary  displacement  sensor, 
bondable  "Krak"  gages,  near-crack  tip  strain  gages,  and  moire  fringe  patterns. 

The  potential  difference  technique  also  was  attempted  on  two  tests,  but  because  of 
unexplained  signals  it  cannot  be  relied  on  at  this  time.  Despite  these  difficulties, 
however,  values  for  Jjc  were  obtained  by  determining  the  threshold  for  crack  initia¬ 
tion.  Table  2  lists  the  specimens  that  were  tested  and  identifies  the  particular 
technique  used  to  obtain  the  information  desired  for  determining  this  toughness 
criterion. 

TEST  PROCEDURE  AND  RESULTS 

A  precrack  had  to  be  placed  in  each  specimen  before  testing.  Each  specimen  was 
installed  in  the  fixture  in  the  usual  manner  and  was  cyclically  loaded  in  the  open¬ 
ing  mode  using  a  hydraulic  actuator-type  machine.  The  applied  load  and  the  corre¬ 
sponding  LPD  of  the  specimen  were  monitored  to  give  the  compliance;  by  combining 
this  with  empirical  relations  in  Reference  25,  the  extent  of  the  crack  growth  was 
determined.  This  operation  was  terminated  when  it  appeared  that  a  suitable  crack 
length  had  been  obtained.  This  fatigue  procedure  was  carried  out  at  two  discrete 
load  ranges.  It  conformed  to  the  procedure  outlined  in  Reference  9,  which  assured 


TABLE  2  -  INSTRUMENTATION  OR  TECHNIQUES  USED  WITH  EACH  SPECIMEN  TO.  OBTAIN  LOAD  POINT 

DISPLACEMENT  AND  DETECTION  OF  CRACK  INITIATION 
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Load  Point  Displacement 
EC  -  Eddy  Current 
FS  -  Fotonic  Sensor 


Crack  Initiation  Detection 

LG  -  Ligament  Gage;  Eddy-Current  Sensor 

NG  -  Notch  Gage;  Eddy-Current  Sensor 

CTC  -  Crack  Tip  Crack  Propagation  Gage 

CTS  -  Crack  Tip  Strain  Gage 

CTK  -  Crack  Tip  "KRAK"  Gage 

MF  -  Mo ire -Fringes 

PD  -  Potential  Drop 


Supplemental  Displacement  Data 
LVDT  -  Linear  Variable 

Differential  Transformer 
HS  -  High  Speed  Film 
S  -  Eddy-Current  Sensor 
in  Specimen  Notch 


®  Number 

1  -  On 

2  -  On 


of  Gages 
One  Side 
Both  Sides 


27 


the  formation  of  a  sharp  crack  tip  that  would  perform  in  a  consistent  manner.  The 


cracking  sequence  finally  selected  initially  called  for  a  5000-lb  cyclical  load 
until  the  crack  reached  approximately  60%  of  the  desired  growth;  then  the  load  was 
reduced  to  2500  lb  to  complete  the  growth  to  about  0.25  in.  (a  =  1.50  in.).  Appro: 
imately  35,000  load  cycles  were  required  to  produce  this  crack.  Figure  12  Is  a 
photograph  of  the  setup  for  this  operation. 


Figure  12  -  Set-Up  in  Cyclic  Load  Machine  for  Fatigue  Cracking 
Compact  Tension  Specimens 


Some  specimens  were  precracked  aC  the  DTNSRDC  Ship  Materials  Engineering 
Department.  Three  discrete  peak  loads  were  used,  starting  with  3500  and  ending 
with  1000  lb  and  upwards  of  100,000  cycles  were  required  to  produce  the  required 
precrack. 

The  fatigue  crack,  rather  than  propagating  at  a  uniform  depth,  tends  to  develop 
in  a  convex  pattern  in  the  direction  of  the  growth  due  to  the  through-thickness 
variation  in  the  stress  state.  This  depth  variation  would  influence  the  test  results, 
and  to  reduce  this  effect,  the  side  grooves,  when  used,  were  machined  after  precrack- 
Ing.  These  side  grooves  were  consistently  0.1  in.  deep  and  had  the  effect  of  reducing 
the  specimen  thickness  to  0.8  in.  at  the  test  section.  Table  3  describes  the  prepa¬ 
rations  made  to  each  specimen  before  testing.  It  identifies  the  final  crack  length 
a  (Including  machined  notch)  from  the  compliance  measurement.  It  also  identifies 
the  specimens  that  contained  side  grooves. 

The  setup  used  in  the  drop  tests  is  shown  in  Figure  13.  In  addition  to  the 
measurements  that  were  needed  to  determine  the  J-integral  for  fracture,  described 
earlier,  the  back  face  strain  (BFS)  was  also  routinely  measured,  mainly  for  use  in 
the  analytical  efforts  (see  Figure  7).  A  block  diagram  of  the  circuitry  used  for 
these  basic  measurements  is  shown  in  Figure  14.  Figure  15  presents  typical  oscillo¬ 
graph  data  obtained  during  the  test.  Peak  values  from  each  channel  of  information 
were  read  routinely;  on  some  channels  the  rise  times  to  peak  values  along  with  other 
pertinent  data  also  were  derived.  Of  course,  to  obtain  the  J-integral  the  applied 
load  and  LPD  Information  needed  additional  analysis,  and  this  was  done  for  most  of 
the  tests  where  the  results  appeared  valid.  These  channels  were  usually  manually 
digitized  and  plots  with  high  res'olution  were  produced  with  use  of  the  Tektronix" 

4081  Graphics  minicomputer.  See  Figures  16  and  17  for  typical  plots.  Through 
recent  software  developments  to  manipulate  arrays  and  perform  integration,  plots  of 


TABLE  3  -  SUMMARY  OF  SPECIMEN  PREPARATION 
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1.47 
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1.56 
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1.53 

NA 

19 

1.57 

1.54 

NA 

5 
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1.53 

NA 

20 
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1.57 

W 

6 

1.47 

1.52 

NA 

21 
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1.48 

NA 

22 
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V 
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13 

1.41 
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V 

29 
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1.51 

W 
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1.42 

1.45 

V 

31 

1.49 

1.51 

W 

32 

1.52 

1.58 

W 

33 

1 .49 

1.49 

W 

'-^See  Figure  9  for  details,  V  -  V  groove,  W  -  Wide  groove, 
NA  -  No  side  grooves. 


Specimens  16,  18-23  were  precracked  by  Ship  Materials  Engi¬ 
neering  Department,  others  by  Structures  Department. 


Figure  16  -  Digitized  Record  of  Applied  Load  versus  Time  for 

Specimen  JC-32 


load  versus  LPD,  and  the  J-lntegral  versus  time  could  also  be  easily  generated. 

(See  Figures  18  and  19  for  examples.)  The  advantage  of  the  plot  shown  in  Figure  19 
is  that  the  critical  value  J  can  be  determined  immediately  If  the  time  for  fracture 
initiation  could  be  identified;  also,  the  time-rate  for  J  (J)  at  fracture  can  be 
approximated  from  the  slope  at  fracture.  [For  justification  of  determining  the  J- 
integral  by  graphical  means,  refer  to  Equation  (5)  in  the  background  section.) 

After  the  test,  each  specimen  was  separated  through  the  ligament  section  to 
reveal  the  crack  faces  so  that  the  specific  fracture  regions  could  be  studied.  To 
help  delineate  these  regions  the  specimen  was  first  heat  tinted;  this  Involved 
placing  it  in  an  oven  and  holding  it  at  a  temperature  of  about  500°F  for  several 
hours.  This  cause/i  oxidation  which  varied  with  surface  roughness.  Thus,  for  fatigue 
fracture,  which  is  characterized  by  a  relatively  smooth  surface,  the  section  would 
appear  lighter  than  the  fractured  section  caused  by  a  single  Impact  as  during  the 
test.  To  separate  the  specimen,  it  was  cooled  to  a  brittle  state  by  holding  it  in  a 
liquid  nitrogen  bath  and  then  placed  over  a  tapered  anvil  that  was  driven  into  the 
front  of  the  notch  by  a  free-fall  Impact  weight.  Figure  20  is  a  photograph  of  the 
fracture  regions  in  a  specimen  that  had  been  severely  tested.  The  zones  are  easily 
discernable.  Each  was  measured  and  the  results  provided  essential  data  for  compu¬ 
ting  the  critical  J-integral.  Then  respective  areas  were  measured  with  a  planimeter 
and,  when  divided  by  the  specimen  thickness,  gave  the  average  lengths  for  each  zone 
(see  Figure  21).  These  measured  precrack  lengths  were  also  tabulated  in  Table  3  for 
comparison  with  those  obtained  from  the  compliance  measurements.  Table  4  is  a  summary 
of  the  magnitudes  for  the  basic  measurements  along  with  the  results  from  analysis  of 
the  data  and  from  measurements  of  the  specimen.  Repeated  tests  were  run  on  some 
specimens  and,  in  those  cases,  only  the  results  from  the  initial  test  are  given. 

Also,  for  specimens  that  were  tested  repeatedly,  the  crack  growth  is  not  meaningful 
and,  therefore,  is  not  listed. 
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versus  Load  Point  Displacement  for  Specimen 


Figure  21  -  Measurements  Typically  Obtained 
After  Separating  Specimen 


FINITE  ELEMENT  ANALYSES 


GENERAL 

In  addition  to  obtaining  sound  experimental  data,  a  major  goal  of  this  program 
was  to  correlate  these  data  with  some  form  of  analytical  solution.  We  originally 
hoped  to  develope  a  reciprocating  relationship  in  which  experimental  data  would  be 
obtained,  and  then  compared  with  data  obtained  by  numerical  analysis.  The  results 
of  this  comparison  would  then  be  used  to  guide  future  testing.  This  approach  would 
be  a  prudent  and  cost  effective  way  of  optimizing  the  experimental  and  analytical 
procedures  at  the  same  time.  However,  the  synergism  of  this  classical  scientific 
approach  did  not  develop  as  anticipated.  The  "short  falls"  Involved  both  the  experi¬ 
ment  (the  problem^  with  data  acquisition  described  in  the  previous  section)  and 
the  newly  evolving  methods  of  analysis.  A  thorough  description  of  these  numerical 
methods,  the  difficulties  encountered  with  their  use,  and  an  evaluation  of  the 
results  follow. 

Several  methods  are  applicable  to  the  analysis  of  a  compact  tension  fracture 
specimen.  In  this  case,  finite  element  computer  codes  were  chosen  as  the  primary 
analytical  method.  With  various  finite  element  computer  codes  it  was  possible  to 
analyze  both  the  dynamic  and  the  static  nonlinear  response  of  the  specimen.  This 
capability  was  crucial  since  the  goals  of  this  program  are  to  determine  the  effects 
of  high  load  rates  on  inelastic  fracture  initiation  and  to  examine  the  applicability 
of  current  nonlinear  static  fracture  criteria  (such  as  the  J-integral)  to  cracked 
specimens  under  dynamic  loading.  The  flexibility  of  general  purpose  finite  element 
programs  was  desirable.  Variations  in  specimen  geometry  (crack  depth),  material, 
and  loading  magnitude  and  rate  could  all  be  easily  modeled  and  altered  for  the  many 
different  test  cases.  Another  factor  influencing  the  choice  of  finite  element  pro¬ 
cedures  as  the  analytical  tool  for  this  program  was  the  large  amount  of  ongoing 
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research  in  the  application  of  finite  elements  to  fracture  problems.  This  research 
would  provide  information  and  Insight  for  these  analyses.  Some  closed  form  theoret¬ 
ical  solutions  were  also  employed  in  conjunction  with  the  finite  elements.  These 
were  Used  primarily  to  verify  and  guide  the  preliminary  static  finite  element 
solutions . 

REVIEW  OF  PRELIMINARY  ANALYSIS 

The  first  analysis  was  performed  on  specimen  JC-6.  The  analysis  of  this  speci¬ 
men  test  has  been  Informally  reported.*  In  the  interest  of  clarity  and  continuity, 
a  summary  of  that  analysis  and  the  insights  gained  are  presented  In  this  section. 

As  experimental  results  indicate,  JC-6  was  among  the  most  severely  tested 
specimens.  It  sustained  a  peak  load  of  7600  lb  (which  was  the  highest  in  the  test 
series)  and  reached  that  peak  load  in  a  very  short  time  (0.6  ms).  With  the  full 
test  series  completed,  it  is  much  easier  now  to  look  back  at  the  results  and  see 
that  JC-6  was  not  a  particularly  good  choice  for  experimental-analytical  comparison. 
As  Table  4  shows,  after  tests  of  JC-7,  it  became  apparent  that  the  tests  with  a 
115-lb  weight  and  36-in.  drop  height  developed  extreme  loadings.  The  weight  and 
drop  height  were  subsequently  lowered  with  the  aim  of  imparting  just  enough  load  to 
achieve  crack  initiation.  The  primary  difficulty  caused  by  the  high  loading  magn¬ 
itude  and  rate  in  JC-6  (and  other  early  series  specimens)  is  decreased  accuracy  and 
reliability  in  the  experimental  results.  The  accurate  ranges  of  the  load  point 
displacement  sensor  and  possibly  other  instrumentation  were  exceeded.  The  limited 
results  available  may  be  reasonably  accurate  (if  the  response  of  the  instrumentation 
was  adequate)  but  the  high  loading  rates  also  cast  doubt  as  to  the  actual  material 


^Rasmussen,  E.  A.,  "Experimental  and  Computational  Analyses  of  a  Dynamically 
Loaded  <"ompact  Tension  Specimen,"  Enclosure  (1)  to  DTNSRDC  letter  lU-1720-59,  dated 
12  May  1981. 
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properties  to  assume  for  the  specimen.  Dynamic  loads  producing  high  rates  of  strain 
can  substantially  affect  yield  stress  as  well  as  other  properties  .26  As  will  be 
seen  later,  the  search  for  an  accurate  description  of  material  properties  is  a  per¬ 
sistent  one. 

One  of  the  most  important  objectives  of  the  comparison  was  to  validate  or  con¬ 
tradict  the  experimentally  obtained  data,  in  which  there  was  already  considerable 
doubt  at  the  time  analysis  was  begun.  For  this  first  analysis,  various  finite  ele¬ 
ment  codes  were  used  to  model  the  specimen.  These  codes  were  ADINA,27  HONDO  11,26 
and  PAPST.16  The  use  of  these  three  codes  allowed  several  different  ideal Izat ion^ 
of  the  specimen  and  thorough  analysis  under  both  static  and  dynamic  loads.  This 
also  served  as  a  means  of  certifying  the  analytical  results  by  comparisons  among 
results  obtained  using  various  codes,  idealizations,  and  types  of  elements. 

The  basic  finite  element  analysis  procedure  follows.  First,  the  compact  tension 
specimen  was  modeled  using  either  4-,  8-,  or  12-noded  isoparametric  elements,  depend¬ 
ing  on  which  code  was  in  use.  The  idealization  for  the  8-noded  element  ADINA  anal¬ 
ysis  is  shown  in  Figure  22.  Two-dimensional  plane  stress  elements  were  employed  in 
this  analysis.  Several  investigators  have  used  finite  element  techniques  to  analyze 
CT  specimens  made  from  various  steels  that  exhibit  elastic-plastic  or  fully  plastic 
fracture.  The  general  outcome  of  these  results  is  that  ( two-din.enslonal )  plane 
stress  elements  yield  displacement  results  that  are  reasonably  close  to  those  ob¬ 
tained  experimentally  or  by  three-dimensional  finite  element  methods .  2*^  » 60  ^o 
special  treatment  such  as  collasped  elements  or  elements  with  an  imbedded  singular¬ 
ity  function  were  used  in  these  cases  to  model  the  crack,  tip  region.  This  type 
of  procedure  was  not  performed  because  the  desired  result  of  the  analysis  was  the 
displacement  at  the  load  point.  Modeling  of  the  stress  singularity  at  the  crack 
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Figure  22  -  ADINA  Two-Dimensional  Finite  Element  Idealization 


tip  was  assumed  to  have  negligible  effect  on  this  displacement.  For  the  static 
condition,  Hilton  and  Glfford^l  confirmed  the  assumption  by  showing  that  various 
treatments  of  the  crack-tip  singularity  yielded  essentially  the  same  J-integral 
value  for  a  finite  element  model  of  an  edge-notched  specimen. 

The  instrumentation  record  in  which  we  had  the  most  confidence  was  the  load 
dynamometer.  For  this  reason,  the  load  was  simulated  by  a  multilinear  curve  and 
used  for  input  in  the  finite  element  solution.  The  load  point  displacement  was  then 
computed  for  the  various  load  steps.  We  hoped  that  these  values  would  be  close  to 
those  measured  by  the  load  point  sensor  in  the  experiment.  From  these  displacement 
values,  load  versus  load  point  displacement  graphs  would  be  developed  so  that  an 
analytical  value  of  the  J-integral  could  be  found  for  comparison  with  the  experiment¬ 
ally  obtained  value.  Also,  from  the  analytical  solution,  stresses  and  plastic  zones 
could  be  determined  that  could  perhaps  give  insight  into  fracture  response  of  speci¬ 
mens  under  the  dynamic  loading. 


Another  aspect  of  the  analytical  solution  was  the  static  analysis  of  the  speci¬ 


men.  We  were  hoping  to  obtain  some  understanding  of  the  dynamic  effects  In  the 
specimen  by  comparing  -a  dynamic  solution  to  a  static  one.  Also,  the  numerical 
solution  could  be  checked  to  determine  Its  convergence  with  the  theoretically  deter¬ 
mined  plastic  limit  load  of  the  compact  tension  specimen.  This  value  was  calculated 
using  the  Merkle-Corten^^  formula  for  the  fully  plastic  load  In  conjunction  with 
the  Green  and  Hundy^^  plastic  constraint  factor  of  1.25.  A  limit  load  value  of 
3788  lb  was  obtained. 

The  most  positive  result  of  this  first  attempt  at  finite  element  analysis  was 
the  close  agreement  between  the  various  programs  used.  As  can  be  seen  in  Figures  23 
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and  24,  under  static  and  dynamic  loading,  load  versus  load  point  displacement  plots 
for  three  different  combinations  of  program  and  element  type  agree  quite  well.  This 
agreement  also  carries  over  to  comparisons  of  the  J-lntegral  and  plastic  zone  sizes 
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Figure  24  -  Specimen  JC-6  Load  versus  Load  Point  Displacement 
for  Dynamic  Finite  Element  Solutions 


as  shown  in  Table  5  and  Figure  25.  This  close  agreement  Instilled  confidence  that 
the  numerical  solutions  wer'  reasonably  correct.  Another  good  agreement,  for  the 
static  loading  condition,  can  be  seen  in  Figure  23.  The  static  solutions  approach 
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the  calculated  limit  load  value  quite  closely;  that  is,  the  deflections  tend  to 
become  extremely  large  as  this  load  value  is  reached.  This  fact  indicates  that  the 
static  finite  element  solutions  are  producing  accurate  displacements. 

As  can  also  be  seen  in  Figures  23  and  24,  the  finite  element  and  experimental 
solutions  do  not  agree  well  at  all.  VJhile  the  accuracy  of  the  finite  element  solu¬ 
tions  was  not  a  known  quantity,  the  repeatability  of  results  by  the  different  codes 
lent  some  credence  to  those  results.  At  the  same  time,  there  was  Inconsistency  in 
some  of  the  experimental  records  (load  point  displacement  being  a  crucial  one). 
These  factors  led  to  a  review  and  modification  of  some  of  the  instrumentation. 
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Figure  25a  -  ADINA  8-Noded  Element  Dynamic  Solution 
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Figure  25b  -  HONDO  Dynamic  Solution 


mm  uii 


III 


!l!!l 


APPLIED  LOADS 


CRACK  TIP-^  1  2  3  4  5 

Figure  25c  -  ADINA  4-Noded  Element  Dynamic  Solution 

Figure  25  -  Dynamic  Finite  Element  Plastic  Zone  Growth 
for  Specimen  JC-6 


SUBSEQUENT  STUDIES 


Further  finite  element  studies  were  performed  after  most  of  the  specimen  tests 
had  been  completed.  At  this  time,  we  had  more  confidence  in  the  Instrumentation, 
which  had  been  fine  tuned  and  debugged  over  the  course  of  the  tests.  Most  important, 
the  load  point  displacement  measurements  obtained  by  the  fiber  optic  sensor  were 
more  consistent  and  reasonable  than  those  obtained  by  eddy  current  sensor.  Thus 
further  finite  element  analysis  was  desirable  at  this  time.  Now  that  the  instrumen¬ 
tation  had  been  refined,  we  hoped  that  further  analyses  would  better  indicate  the 
ability  of  these  methods  to  predict  specimen  response  under  dynamic  loading. 

Using  the  same  general  procedure  used  in  the  preliminary  analyses,  four  more 
specimens  were  analyzed.  The  specimens  had  markedly  different  crack  depths  and  load 
histories.  Figure  26  shows  the  specimen  numbers  and  crack  depths  along  with  the 
respective  load  histories  of  the  specimens  used  in  this  second  set  of  finite  element 
analyses.  The  load  histories  shown  are  the  actual  input  data  for  the  computer  solu¬ 
tions.  For  this  reason,  the  curves  are  represented  in  a  piece-wise  linear  fashion. 
These  Idealized  load  curves  were  taken  directly  from  the  experimental  load  records 
but  with  an  effort  to  generally  smooth  the  curves  and  disregard  extraneous  spikes. 

The  slower,  lower  magnitude  loading  had  fewer  spikes  and  required  less  smoothing. 

The  load  history  for  JC-6,  analyzed  previously,  is  also  provided  in  Figure  26  for 
comparison.  As  can  be  seen,  the  histories  of  the  four  later  tests  exhibit  loads  of 
lower  magnitude  and  much  longer  duration.  This  change  in  load-history  character  is 
due  to  the  lighter  drop  weight,  lower  drop  height  and,  for  specimens  JC-23  and  JC-27, 
a  different  load  programmer.  The  slower  rate  of  these  curves  made  the  finite  element 
analysis  somewhat  easier  to  perform,  particularly  the  choice  of  time  steps  (increment 
between  solutions)  and  the  providing  of  more  complete  solutions. 


The  ADINA  program  and  the  Idealization  shown  In  Figure  22  were  used  exclusively 
for  this  series  of  analyses.  This  combination  had  yielded  the  best  results  In  the 
previous  analysis  of  specimen  JC-6  dynamic  loading. 

There  were  some  modeling  difficulties  In  the  second  series  of  analyses  due  to 
the  side  grooves  on  the  actual  specimens;  the  grooves  reduced  the  thickness  from  1 
to  0.8  in.  In  the  area  of  the  crack  and  uncracked  ligament.  These  side  grooves  were 
of  two  types,  either  1/8-ln.  radius  curves  or  V-notches.  Modeling  this  variable 
thickness  In  the  specimen  (without  using  an  expensive  three-dimensional  dynamic 
solution)  required  some  manipulation  of  the  previously  used  finite  element  model. 

The  side  grooves  were  modeled  by  reducing  the  thickness  of  the  bottom  row  of  elements 
in  the  Idealization  (Figure  22)  from  1  to  0.8  in.  This  change  was  Important  because, 
for  two-dimensional  plane  stress  elements,  the  element  thickness  Is  a  component  of 
the  stiffness.  The  actual  and  idealized  specimens  are  shown  In  Figure  27.  For  speci¬ 
mens  with  wide  grooves,  about  the  same  volume  was  removed  In  the  side  grooved  region 
in  both  cases.  Maintaining  this  equivalence  of  volume  was  considered  important  so 
that  about  the  same  mass  would  be  present  in  the  region  of  the  crack  tip. 

The  load  pin  hole  was  not  modeled.  Unlike  a  static  compact  tension  test,  which 
uses  a  loosely  fitting  pin  to  apply  load,  the  hole  was  almost  completely  filled  with 
a  tightly  fitted  bolt.  Leaving  the  Idealized  specimen  without  a  pin  hole  models 
the  mass  effect  of  this  bolt-specimen  interface  reasonably  well.  Consideration  of 
the  load  pin  and  hole  led  to  another,  perhaps  more  pertinent  question:  Is  it  neces¬ 
sary  to  model  the  stiffness  and  mass  of  the  bottom  loading  clevis,  since  the  applied 
load  is  measured  at  the  end  of  the  clevis  opposite  the  specimen?  There  was  a  consis¬ 
tent  lag  in  displacement  response  behind  load  response  in  the  experimental  records. 

The  exact  cause  of  this  time  lag  was  not  clear,  but  we  decided  that  the  experimental 
displacement  response  curves,  shifted  to  eliminate  this  lag,  were  the  most  appropriate 
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Figure  27  -  Actual  and  Idealized  Side  Groove  Geometry 

goals  for  the  finite  element  analyses.  If  agreement  with  these  was  achieved,  then 
the  relative  effect  of  the  clevis  could  be  further  quantified.  This  will  be 
addressed  later  in  this  section. 

A  method  of  verification  was  needed  to  determine  the  validity  of  the  side 
grooved  specimen  model.  A  comparison  of  the  static,  linear  elastic  finite  element 
solution  to  that  obtained  by  a  compliance  equation  would  provide  this  verification 
In  the  elastic  range.  Compliance  is  defined  in  the  linear  elastic  case  as  the  load 
point  displacement  per  unit  load.  The  equation  chosen  was  developed  by  Holsberg.^S 
It  relates  the  compliance  to  the  crack  depth  (a)  and  total  specimen  width  (W): 


where  C  =  Compliance  =  load  point  displacement/load 

y  =  Poisson's  ratio  =0.3 

E  =  Young's  modulus  =  29.0  ^  10^  psl 

B  =  Specimen  thickness  =  0.8  In.  (with  side  grooves) 

1.0  in.  (without  side  grooves) 

a  =  Crack  depth  =  1.51  In. 

W  =  2.0  In. 

For  the  side  grooved  specimen  (B  =  0.8  in.)  the  equation  yields; 

C  =  7.791  X  10-6  In. /lb 

and  for  a  similar  specimen  without  side  grooves  (B  =  1.0  In.); 

C  =  6.232  X  10-6  In. /lb 

Two  finite  element  analyses  were  performed:  one  on  a  specimen  with  the  side 
grooves  and  one  without  side  grooves.  For  linear  elastic  behavior,  the  results  are 
shown  compared  to  those  obtained  by  the  compliance  equation  in  Table  6.  As  can  be 

TABLE  b  -  LOAD  POINT  DISPLACEMENTS  BY  COMPLIANCE  ATJD  FINITE 
ELEMENT  METHODS  FOR  SIDE  GROOVED  AND 
FULL  THICKNESS  SPECIMENS 


Ef  feet  ive 

Spec imen 

Width 
(  in . ) 

Load  Point  Displacement  Under 

Static  lOOO-lb  Load 
(in.) 

Finite  Element 

Compliance  Equation 

0.8 

(Side  Groove) 


0.00710 


0.00779 


seen,  the  results  for  the  1 .0-in. -thick  specimen  agree  quite  well  with  the  finite 
element  solution.  The  results  for  the  0.8-in. -thick  side  grooved  solution  show  that 
the  finite  element  solution  is  considerably  stlffer  (10%)  than  that  obtained  by  the 
compliance  equation.  It  is  not  clear  whether  this  discrepancy  is  due  to  inaccura¬ 
cies  in  the  finite  element  model  or  the  assumption  of  0.8  in.  for  the  specimen  thick¬ 
ness  in  the  equation  producing  an  overestimation  of  specimen  compliance.  It  was  felt 
that  for  the  wide  side  grooved  specimen,  the  compliance  should  be  close  to  that  of  a 
0.8  in.  constant  thickness  specimen.  Obtaining  this  agreement  would  require  a  re¬ 
duction  in  stiffness  of  the  0 .8-ln. -thick  elements. 

The  most  obvious  way  to  reduce  the  stiffness  is  to  further  decrease  the  thick¬ 
ness  of  the  elements  to  less  than  0.8  in.  This  was  undesirable  because  it  would 
produce  a  modeled  crack  width  smaller  than  that  of  the  actual  specimen.  This  would 
lead  to  an  overprediction  of  displacements  for  loads  producing  plasticity  in  the 
crack  region.  A  better  possibility  was  to  reduce  the  Young's  modulus  (E)  of  the 
thinner  elements.  The  only  problem  caused  by  changing  the  Young's  modulus  would  be 
a  variation  iji  the  amount  of  elastic  strain  occurring  before  the  onset  of  plasticity. 
If  the  new,  smaller  Young's  modulus  is  close  to  the  old  value,  there  would  be  a  very 
minor  difference  with  plasticity  occurring  at  perhaps  0.1%  higher  elastic  strain. 
With  displacements  of  the  specimen  being  the  primary  results,  and  considering  finite 
element  accuracy,  this  slight  error  was  thought  to  be  insignificant. 

Having  decided  that  reducing  the  Young's  modulus  of  the  thinner  element  was  the 
preferred  method  to  match  the  compliance  solution,  we  had  to  determine  what  reduced 
value  was  best.  Using  the  same  linear-elastic  analysis  as  before,  various  values 
for  Young's  modulus  were  substituted  into  the  finite  element  solution.  The 
results  are  shown  in  Table  7  and  are  compared  to  the  compliance  equation  results. 


TABLE  7  -  LOAD  POINT  DISPLACEMENT  BY  COMPLIANCE  AND 
FINITE  ELEMENT  METHODS  FOR  A  SIDE 
GROOVED  SPECIMEN  VARYING 
YOUNG'S  MODULUS 


Load  Point  Displacement  Under  Static 

1000-lb  Load  (in.) 

Compliance 
Equat ion 

E=29  x  10  psi 

Finite  Element  Method 

E=27  X  10^  psi 

E=26.1  X  10^  psi 

E=24  X  10^  psi 

0.00779 

0.00732 

0.00746 

0.00777 

An  elastic  modulus  of  24  ^  10^  psl  produced  load  point  displacements  closely 
matching  those  predicted  by  the  compliance  method. 

After  developing  a  suitable  model  for  the  side-grooved  compact  tension  specimen, 
the  next  step  was  to  conduct  dynamic  analyses  using  the  experimental  load  versus 
time  histories  as  input.  Specimen  JC-27  was  analyzed  first.  As  can  be  seen  in 
Figure  26,  the  load  history  for  JC-27  was  smoothly  shaped  and  of  relatively  long 
duration.  Peak  magnitude  was  4150  lb  at  a  time  of  approximately  2.0  ms  (these  may 
differ  slightly  from  actual  experimental  values  due  to  piecewise  linear  approxima¬ 
tion  of  load  histories). 

Choosing  accurate  material  properties  for  the  dynamic  analyses  caused  some  diffi¬ 
culty.  Only  one  static  coupon  test  had  been  performed  on  the  HY-80  plate  from  which 
all  the  specimens  were  machined.  This  test  produced  a  0.2%  offset  yield  stress  of 
88,380  psl  and  an  engineering  ultimate  stress  of  105,300  psl  at  approximately  6% 
strain.  The  ADINA  finite  element  program  is  designed  so  that,  for  an  elastic- 
plastic,  isotropic  hardening  material,  only  a  bilinear  stress  strain  curve  may  be 
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specified.  The  required  input  parameters  are  Young's  modulus,  Poisson's  ratio, 
yield  stress,  and  strain  hardening  modulus.  These  requirements  limit  the  accurate 
modeling  of  the  nonlinear  material.  Choosing  an  accurate  strain  hardening  modulus 
to  best  represent  the  nonlinear  behavior  of  the  material  is  particularly  difficult. 
Hammel  et  al.29  addressed  this  problem  for  a  statically  loaded  steel  compact  tension 
specimen,  and  their  best  results  were  obtained  using  a  strain  hardening  value  tangent 

to  the  early,  steeper  portion  of  the  hardening  curve. 

After  considering  that  work  and  the  error  induced  by  having  only  one  material 
property  test,  we  decided  to  input  the  88,380  psi  value  as  yield  stress  and  to  use  a 
hardening  modulus  of  300,000  psi,  which  is  indicative  of  the  slope  of  the  true 
stress/true  strain  curve  between  the  values  obtained  from  the  coupon  test.  The 
Young's  modulus  (for  other  than  the  side-grooved  region)  and  Poisson's  ratio  were 
29.0  X  10^  psi  and  0.3  respectively.  The  idealized  stress-strain  curve  is  shown 
in  Figure  28. 

Specimen  JC-27  was  analyzed  under  plane  stress  conditions.  The  primary  desired 
outcome  from  this  and  subsequent  analyses  was  agreement  with  the  experimental  load 
point  displacement  response.  If  good  agreement  occurred  between  these  displacements, 

other  calculated  values,  such  as  the  J-integral  would  also  be  in  close  agreement 

because  the  experimental  load  history  is  used  as  input  for  the  finite  element  anal¬ 
yses.  For  this  first  case,  the  analytically  obtained  displacments  were  much  larger 
than  those  from  the  experimental  record.  It  was  determined  that  these  large  displace¬ 
ments  could  be  at  least  partly  due  to  the  use  of  static  material  properties  in  the 
dynamic  solution.  The  decision  was  made  to  modify  the  material  properties  in  order 
to  obtain  a  solution  in  agreement  with  the  experimental  results.  Yield  stress  was 
the  most  obvious  choice  because  there  are  a  limited  number  of  references  available 
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documenting  measured  Increases  in  HY-80  yield  stress  with  Increasing  strain 
rate. 26. 33  -phe  effect  of  strain  rate  on  other  material  properties  of  HY-80  could 
not  be  found  in  the  available  literature. 

Various  Increases  in  yield  stress  were  then  used  in  an  effort  to  decrease  the 
deformation  of  the  Idealized  specimen.  As  shown  In  Figure  29,  calculated  displace¬ 
ments  were  nearly  inversely  proportional  to  yield  stress.  A  value  of  110,480  psi  (a 
25%  Increase  over  88,380  psi)  produced  displacement  results  very  close  to  those 
obtained  from  the  experiment.  This  yield  stress  value  seemed  quite  large  even  with 
the  relatively  high  strain  rates  present.  The  measured  back  face  strain  rate  for 
this  test  was  on  the  order  of  12  s“l.  According  to  References  26  and  33,  tests 
conducted  on  HY-80  at  approximately  these  strain  rates  Indicated  the  yield  stress 
Increase  would  be  more  on  the  order  of  5  to  10%.  The  back  face  strain  rate  Is  not 
representative  of  the  strain  rate  near  the  crack  tip,  but  much  of  the  plasticity 
present  in  the  specimen  does  occur  at  the  back  face  so  the  strain  rate  measured 
there  is  important  in  terms  of  contribution  to  overall  specimen  deformation.  It 
seems  likely  that  it  is  not  only  the  material  properties,  but  also  the  finite 
element  model  contributing  to  the  overly  flexible  analytical  response.  In  an 
effort  to  resolve  this  material  property  question,  more  analyses  were  performed 
using  other  experimental  records. 

The  JC-23  specimen  test  was  particularly  suitable  for  finite  element  analysis 
because  the  load  was  of  relatively  small  magnitude  and  was  applied  over  quite  a  long 
duration  (Figure  26).  An  extensive  analysis  was  performed  on  this  specimen  in  hopes 
of  resolving  the  yield  stress  question.  In  addition,  we  investigated  other  methods 
of  modeling  the  geometry  and  loading  of  the  specimen  to  better  simulate  the  experi¬ 
mental  case.  The  first  analyses  were  done  using  a  constant  hardening  modulus  of 
300,000  psi.  The  yield  stress  was  varied  to  obtain  displacement  results  matching 
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Figure  29  -  Experimental  and  Finite  Element  Load  Point  Displacement 
versus  Time  Results  for  Specimen  JC-27 
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those  from  the  experimental  record  (the  procedure  was  similar  to  that  used  In  the 
JC-27  analysis).  The  results  of  this  Investigation  are  shown  in  Figure  30.  The 
displacements  for  a  yield  stress  of  92,800  psl,  a  5%  increase  over  the  static  meas¬ 
ured  yield,  are  quite  close  to  the  experimental  deflection  curve.  This  figure  is 
also  a  good  example  of  the  sensitivity  of  the  finite  element  displacements  to  yield 
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Figure  30  -  Experimental  and  Finite  Element  Load  Point  Displacement  versus 

Time  Results  for  Specimen  JC-23 


stress  change.  As  shown,  load  point  displacement  differed  as  much  as  40%  with  a 
yield  stress  change  of  only  5%  (88,380  to  92,800  psl).  This  Important  result  Indi¬ 
cates  that  accurate  measurement  of  the  material’s  yield  stress  Is  crucial  for  this 
type  of  testing.  In  this  program,  for  which  only  one  coupon  test  was  performed,  the 
measured  yield  stress  Is  surely  accurate  to  within  5%  at  best.  These  results  show 
clearly  that  this  accuracy  Is  Inadequate.  As  with  JC-27,  JC-23  undoubtedly  experi¬ 
enced  some  strain  rate  effects.  Figure  26  shows,  however,  that  the  loading  rate  for 
JC-23  was  less  than  that  for  JC-27.  Also,  from  the  experimental  records,  the  back- 
face  strain  rate  was  8.0  which  is  less  than  that  measured  for  JC-27.  Yield 

stress  Increases  for  HY-80  under  strain  rates  of  this  magnitude  should  be  only  on 
the  order  of  3  to  5%. 28, 33  Thus  for  this  test  (and  probably  for  the  lower  load 
rates  in  general),  pinpointing  the  amount  of  yield  stress  increase  Is  not  as  impor¬ 
tant  as  obtaining  an  accurate  static  yield  stress  for  the  material. 

Perhaps  the  most  Important  result  of  these  analyses  thus  far  is  that  once  the 
correction  of  material  properties  is  performed,  the  finite  element  results  model  the 
experimentally  obtained  displacements  quite  well.  A  predictive  capability,  there¬ 
fore,  is  certainly  obtainable.  Also,  the  Idealization  chosen  Is  probably  adequate 
for  this  analysis. 

To  further  verify  these  results,  an  effort  was  made  to  fine-tune  the  idealiza¬ 
tion  by  modeling  the  contribution  of  the  lower  loading  clevis.  As  mentioned  earlier, 
the  effect  of  the  stiffness  and  mass  of  the  clevis  on  the  specimen  was  unknown. 

Since  the  applied  load  is  measured  at  the  necked-down  area  of  the  lower  clevis,  the 
clevis  should  be  modeled  as  an  extension  of  the  specimen.  The  mass  and  stiffness 
were  modeled  by  means  of  a  linear  elastic  truss  element  and  a  concentrated  nodal 
mass.  The  neck  of  the  clevis,  the  flexibility  of  which  was  thought  to  be  a  possible 
cause  of  the  consistent  lag  of  experimental  displacement  measurements  behind  the  load 
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record,  was  modeled  by  a  truss  element  having  the  same  length  and  stiffness.  The 
rest  of  the  clevis  was  assumed  to  be  rigid  and  was  treated  as  a  concentrated  mass  at 
the  specimen  load  point  as  shown  In  Figure  31.  Most  likely,  this  model  would  accu¬ 
rately  reproduce  any  deformat ional  and  inertial  effects  of  the  actual  clevis  on  the 
measured  load-point  displacement.  This  idealization  was  run  using  the  same  JC-23 
input  values  that  produced  the  best  results  shown  in  Figure  30.  The  load  point 
displacement  obtained  from  this  analysis  is  shown  in  Figure  32;  the  experimental 


ACTUAL  CLEVIS  GEOMETRY 


IDEALIZED  MODEL  OF  CLEVIS  INCLUDED  IN  IDEALIZATION  OF  SPECIMEN 

Figure  31  -  Finite  Element  Model  of  Loading  Clevis 
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value  and  the  best  result  without  modeling  the  clevis  are  both  included.  As  shown 
in  the  figure,  the  clevis  model  had  a  slight  damping  effect  and  reduced  the  displace¬ 
ment  a  small  amount.  No  significant  time  lag  was  evident,  but  the  first  computer 
displacement  output  was  probably  too  late  for  this  to  be  a  conclusive  result.  A 
second  attempt  was  made,  this  time  including  the  approximate  stiffness  of  the  body 
of  the  clevis  (not  just  the  neck)  in  the  truss  element.  This  more  detailed  clevis 
model  was  used  in  conjunction  with  the  JC-27  idealization  and  loading  input. 

Specimen  JC-27  exhibited  a  displacement  lag  of  about  0.3  ms  behind  the  load  record 
when  tested.  We  hoped  that  the  more  accurate  model  of  the  clevis  and  an  additional 
finite  element  displacement  output  at  0.2  ms  would  show  evidence  of  this  lag.  The 
results  of  the  JC-27  analysis  are  presented  in  Figure  33.  Shown  are  the  displacement 
history  obtained  without  modeling  the  clevis  and  the  displacements  for  the  case  with 
the  accurately  modeled  clevis.  Close  examination  of  the  results  for  a  time  of  0.2 
ms  reveals  that  only  a  very  small  displacement  (perhaps  not  dlscernable  on  the  exper¬ 
imental  records),  had  occured  up  to  this  time.  This  delay  caused  by  the  clevis 
flexibility  could  be  a  major  component  of  the  experimental  time  lag.  In  this  case, 
the  analysis  did  not  clearly  show  the  full  0.3  ms  delay  of  the  experimental  displace¬ 
ments,  but  it  must  also  be  remembered  that  no  flexibility  in  the  load  bolt  connection 
was  Included. 

The  fact  that  the  experimental  time  lag  appears  to  be  largely  due  to  the  flexi¬ 
bility  of  the  loading  clevis  (more  correctly,  due  to  the  load  measurement  taking 
place  on  the  clevis  and  not  at  the  load  pin)  is  an  Important  finding.  It  means 
there  were  no  large  losses  in  the  system,  only  a  delay  of  the  force  acting  on  the 
specimen  Itself.  Although  the  mass  of  the  clevis  does  appear  to  play  a  slight  role 
in  determining  magnitude  of  load  point  displacement  in  the  two  analyses  performed 
(Figures  32  and  33),  there  does  not  appear  to  be  a  specific  trend.  Considering  the 
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number  of  other  uncertainties  and  observing  that  the  inertial  effect  of  this  mass 
tends  to  be  small,  it  appears  reasonable  to  Ignore  the  inertial  contribution  of  the 
clevis . 

■These  findings  justify  the  shifting  of  the  experimental  displacement  curves  so 
that  displacement  coincides  with  load  increase.  Unless  the  experimental  displacement 
curve  is  shifted,  the  results  obtained  for  the  load  vercus  load  point  displacement 
curve  and,  subsequently,  the  J-integral  will  be  incorrect.  The  specimen  will  appear 
considerably  stlffer  than  it  actually  is  because  of  the  apparent  lack,  of  displacement 
during  the  early  time  loading.  The  finite  element  results  indicate  that  this  lack 
of  displacement  is  caused  primarily  by  the  flexibility  of  the  clevis  so  that  it  is 
reasonable  to  remove  the  clevis,  in  effect,  and  apply  the  load  directly  to  the  speci¬ 
men.  Shifting  the  experimental  displacement  curve  suggests  that  there  is  no  clevis 
and  that  the  only  response  is  the  response  of  the  specimen.  These  conclusions  indi¬ 
cate  that  the  best  method  for  reproducing  shifted  displacement  response  of  the  spec¬ 
imen  is  applying  the  load  directly  to  the  idealized  specimen  and  not  modeling  the 
clevis  at  all.  This,  of  course,  is  true  only  if  the  inertial  effects  of  the  clevis 
are  small  (as  they  appear  to  be).  Thus,  we  apparently  had  completed  a  full  circle. 
The  first  analyses  were  performed  on  JC-27  and  JC-23  without  modeling  of  the  load 
clevis.  It  now  appears  that  this  type  of  analysis  is  most  appropriate  for  reproduc¬ 
ing  the  shifted  experimental  displacement  curves. 

Having  determined  the  most  appropriate  idealization  geometry,  we  decided  to 
continue  finite  element  specimen  analysis  to  further  investigate  the  material  prop¬ 
erty  questions  raised  earlier.  Two  more  specimen  tests  were  Investigated.  Specimens 

\ 

JC-15  and  JC-22  were  chosen  because  they  were  subjected  to  loading  that  was  basically 
different  from  the  loading  of  the  previously  analyzed  specimens.  Their  load 
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histories  are  shown  In  Figure  26  along  with  the  load  histories  of  the  other  analyzed 
specimen  tests.  The  loading  rates  for  JC-15  and  JC-22  were  somewhat  higher  at  the 
beginning  of  the  load  history  curve,  then  peaked  quickly,  and  the  load  value  oscil¬ 
lated  once  before  declining  again.  This  shape  resulted  from  the  properties  of  the 
load  programmer  used  for  these  tests — a  nylon  plug,  which  produced  load  characteris¬ 
tics  considerably  different  from  those  of  the  1/8-ln. thick  lead  pad  used  In  tests 
of  the  previously  analyzed  specimens.  Specimens  JC-15  and  JC-22  also  underwent  a 
slight  amount  of  crack  extension,  indicating  that  a  critical  value  of  the  J-lntegral 
(die)  reached  in  these  tests,  although  hopefully  not  exceeded  by  too  large  an 
amount.  We  wanted  to  analyze  specimens  in  which  crack  growth  had  occurred  to  find 
any  notlcable  differences  between  the  experimental  displacement  curve  and  the  finite 
element  results  due  to  the  crack  extension  that  would  not  be  accounted  for  in  the 
analysis . 

First,  JC-22  was  analyzed  using  the  same  basic  procedure  described  earlier. 
Because  of  the  load  rate  and  magnitude,  the  JC-22  specimen  behavior  ought  to  be 
similar  to  that  of  the  JC-23  specimen.  As  with  the  previous  analyses,  material 
yield  stress  had  to  be  varied  to  produce  finite  element  load  point  displacements 
approximating  those  measured  In  the  experiment.  As  shown  In  Figure  34,  a  yield 
stress  of  98,990  psl  produced  finite  element  displacements  reasonably  close  to  the 
experimental  results.  This  value  was  12%  above  the  static  yield  as  compared  to  the 
5  and  25%,  respectively,  required  to  achieve  reasonably  close  solutions  In  specimens 
JC-23  and"  JC-27.  Because  the  back  face  strain  rate  for  this  test  was  not  obtained, 
no  strain  rate  comparison  was  available.  Figure  26  shows  that  the  Initial  load  rate 
Is  somewhat  higher  than  ior  JC-23  and  about  the  same  as  for  JC-27.  From  the  results 
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of  the  three  analyses  performed  thus  far,  it  appears  that  the  smaller  experimental 
displacement,  accounted  for  in  the  analysis  by  increased  yield  stress,  is  load  rate 
(strain  rate)  dependent. 

A  pleasing  product  of  the  JC-22  solution  was  the  manner  in  which  the  finite 
element  solution  followed  the  oscillation  of  the  experimental  displacement  curve. 
This  agreement  indicates  that  the  effect  of  the  lower  Young's  modulus  in  the  side 
groove  region  seems  to  be  minimal,  since  the  displacement  during  the  unloading  and 
reloading  agrees  with  the  experimental  results  quite  well.  The  nagging  question  of 
why  the  large  Increase  in  yield  stress  is  necessary  to  obtain  an  analytical  solution 
approximating  the  experiment  remains.  As  was  discussed  previously,  strain  rate 
effects  are  a  possible  explanation  especially  as  they  may  affect  the  hardening  modu¬ 
lus.  In  the  case  of  JC-22,  if  the  strain  rates  are  on  the  order  of  those  found  in 
JC-27,  this  could  explain  nearly  all  of  the  required  yield  stress  increase.  However 
considerable  doubt  remains. 

Specimen  JC-15  was  the  most  severely  tested  specimen  in  this  set  of  four  anal¬ 
yses.  It  was  also  the  only  specimen  of  the  four  having  narrow  V-shaped  side  grooves 
As  shown  in  Figure  26,  the  initial  load  rate  and  peak  load  were  quite  high.  From 
the  previous  cases,  one  would  expect  at  least  a  25%  increase  in  yield  stress  to  be 
required  to  obtain  displacements  near  those  measured  in  the  experiment.  The  plots 
of  load  point  displacement  versus  time  obtained  from  the  experiment  and  from  the  ana 
lytlcal  solutions  are  compared  in  Figure  35.  In  this  case,  a  40%  yield  stress  in¬ 
crease  was  needed  to  obtain  displacements  relatively  near  the  experimental  values. 
The  back  face  strain  record  for  this  test  indicated  a  strain  rate  of  approximately 
20  s"l.  This  strain  rate  could  cause  up  to  a  10%  Increase  in  the  yield  stress  for 
HY-8033  but  certainly  not  the  40  or  even  50%  Increase  required  to  match  the  experi¬ 
mental  displacements. 
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Figure  35  -  Experimental  and  Finite  Element  Load  Point  Displacement  versus 

Time  Results  for  Specimen  JC-15 
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The  finite  element  solution  for  JC-15  Indicates  that  the  Idealization  used, 
specifically  the  method  In  which  the  aide  grooves  were  modeled,  was  Inappropriate 
for  this  analysis.  Intuitively,  one  would  expect  the  specimen  with  the  V-shaped 
grooves  to  be  somewhat  less  compliant  elastically  than  the  finite  element  model 
used.  Apparently,  It  Is  also  considerably  less  plastically  compliant  Chan  Che  finite 
element  Idealization.  The  finite  element  model  overpredlcCs  the  actual  displacements 
because  the  reduced  thickness  section  of  the  model,  which  was  developed  to  represent 
the  wide  side  grooves,  does  not  accurately  represent  the  V-shaped  grooves  present  in 
specimen  JC-15  (see  Figures  9  and  27).  In  effect,  the  plastic  zones  at  the  crack 
tip  and  back  face  in  the  reduced  thickness  region  of  the  model  expand  further  under 
a  given  load  Chan  Is  appropriate.  This  Is  also  a  possible  explanation  for  the 
overprediction  of  displacements  that  occurred  In  the  analysis  of  specimen  JC-27.  As 
shown  in  Table  8,  which  gives  pertinent  experimental  and  finite  element  data  for  the 
four  specimens  analyzed,  JC-27  sustained  the  highest  load  of  the  wide  side  grooved 
specimens.  At  higher  loads,  possibly  producing  larger  plastic  zones,  the  side  groove 
model  becomes  Increasingly  Inaccurate  because  the  gradual  increase  In  thickness  of 
the  actual  specimen  Is  not  modeled.  This  means  that  the  plastic  zones  will  be  larger 
In  the  finite  element  analysis  than  In  the  actual  specimen  causing  an  overprediction 
of  displacements.  The  smaller  magnitude  of  the  loads  on  specimens  JC-22  and  JC-23 
may  have  reduced  this  effect  to  the  point  that  accurate  answers  could  be  obtained 
using  reasonable  dynamic  yield  stresses. 

DISCUSSION 

Over  the  course  of  ,the  experimental  program,  an  accurate  method  for  testing 
compact  tension  specimens  under  dynamic  loading  was  developed.  This  capability 
raises  questions  of  how  the  data  obtained  by  this  method  should  be  applied,  As  a 
means  of  comparison,  the  J-lntegral  was  routinely  determined  from  the  specimen  test 
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load  and  load  point  displacement  data.  These  values  were  then  used  as  a  measure  of 
the  material  Jjc  by  determining  which  calculated  value  of  J  best  corresponded  to 
crack  Initiation.  One  must  question  the  application  of  standard  static-test  methods 
of  determining  the  J-lntegral  In  light  of  (1)  the  Inertial  effects  caused  by  the 
dynamic  loading  and  (2)  the  associated  vibratory  loading  and  unloading  of  the  speci¬ 
men.  The  actual  significance  of  these  effects  cannot  easily  be  quantified.  It  Is 
possible  that,  for  these  specimens,  Inertial  effects  are  not  that  large,  considering 
the  damping  caused  by  the  programmers  and  load  line  velocities  of  10  to  50  In./s. 

The  relatively  smooth  load  histories  obtained  In  the  test  data  seem  to  support  this 
contention.  Also,  the  J-lntegral  was  the  only  readily  available  measure  of  material 
toughness  that  could  be  easily  applied  to  these  tests.  Thus  the  J  and  Jjc  values 
were  used  for  comparisons,  although  doubt  about  their  ultimate  validity  has  not  been 
resolved. 

The  J-lntegrals  obtained  for  specimen  tests  where  the  Instrumentation  seemed 
to  be  accurate  are  presented  In  Table  9.  This  table  Is  separated  Into  two  sections: 
specimens  tested  at  load  rates  produced  using  a  nylon  load  programmer  and  those 
tested  using  a  lead  load  programmer.  The  J-lntegrals  were  calculated  by  determining 
the  maximum  area  under  the  load  versus  LPD  curve  and  applying  the  Merkle-Corten 
equation. No  provision  was  made  to  account  for  crack  extension,  so  that  the 
calculated  Jmax  be  Incorrect  for  specimens  having  undergone  substantial  crack 
growth.  Considering  the  Instrumentation  problems  incurred  in  this  program,  it  Is 
especially  pleasing  to  observe  how  well  these  test  results  appear  to  bracket  a  crit- 

O 

leal  Jj  value  of  approximately  1000  in.-lb/ln.  for  both  loading  situations.  This 
Is  slightly  lower  than  the  1110  In.-lb/in.*  value  obtained  from  three  static  J^^ 
tests  on  the  same  plate,  although  with  the  relatively  few  available  data  points. 


TABLE  9  -  SUMMARY  OF  EXPERIMENTAL  Jic  RESULTS 


Specimen 

Drop  Height 
(in.) 

Peak  Load 
(lb) 

Max  Ji 

(in. -lb/ in. ^) 

Crack  Extension 
(in.) 

Nylon  Programmer 

JC-10 

6 

4410 

7  60 

None 

JC-11 

5 

4800 

720 

None 

JC-14 

7 

5325 

1070 

None 

JC-15 

8 

4690 

1040 

<0.01 

1 

JC-20 

8 

2940 

1000 

■k 

JC-22 

8 

2950 

1120 

0.02 

1 

Lead  Programmer 

JC-2  3 

12 

2400 

670 

None 

JC-2  6 

12 

4100 

1500 

0.02 

JC-27 

10 

4080 

960 

None 

JC-29 

10 

4200 

1040 

0.01 

JC-32 

12 

3250 

1650 

0.04 

JC-33 

11 

1 

3100 

_ 1 

1380  ' 

1 

0.02 

*Speclmen  was  tested  twice;  visible  extension  observed  after  it  was  broken 
open . 


the  difference  is  probably  not  significant.  Hasson  and  Joyce^2  reported  a  similar 
result  for  compact  specimens  loaded  in  a  servohydraullc  test  machine.  In  their 
work,  for  HY-80  steel,  as  determined  by  the  Merkle-Corten  formula,  was  comparable 

(at  room  temperature)  for  load  point  rates  of  5.1  ln,/s  and  1 .6  10“^  In./s.  At 

lower  than  room  temperature,  the  rapidly  loaded  specimens  produced  Jjq  values  lower 
than  those  found  under  static  conditions. 

The  experimental  results  obtained  after  much  alteration  and  fine  tuning  of 
instrumentation  seem  to  be  reasonably  believable  and  accurate.  Questions  remain  as 
to  what  fracture  toughness  assessment  is  most  applicable  in  this  case  considering 
the  apparent  shortcomings  of  the  J-integral.  New  toughness  parameters,  such  as  J’ 
and  AT,  which  are  valid  in  the  dynamic  and  nonproportional  loading  cases,  are 
being  investigated.^^ >35  Though  not  readily  available  at  this  writing,  these 
or  other  parameters  may  allow  an  ultimate  assessment  of  this  work. 

The  finite  element  analyses,  we  had  hoped,  would  play  a  role  in  the  character¬ 
ization  of  the  material  toughness.  The  problems  encountered  in  matching  the  load 
point  displacement  histories  of  the  experiment  diminished  that  hope.  The  large 
yield  stresses  required  to  retard  the  finite  element  load  point  displacements  require 
more  investigation.  There  are  several  possible  explanations  for  the  discrepancies 
between  the  finite  element  results.  The  first  explanation  presumes  inaccuracies 
remain  in  the  load  point  displacement  measurement.  This  appears  unlikely  in  tests 
where  the  fiber  optic  Fotonic  sensor  was  used  for  displacement  measurement.  Much 
"double  checking"  was  done  with  this  sensor  including  measurement  of  set  displace¬ 
ments  with  the  sensor  and  feeler  gages.  Results  of  this  verification  Indicated  that 
•  \ 

the  load  point  sensor  is  probably  accurate  to  within  10%. 


A  second  possible  explanation  Is  that  the  finite  element  solution  was  dependent 
on  a  poor  model  of  the  material  properties,  A  bilinear  stress-strain  model  Is  not 
as  accurate  as  power  hardening  or  multilinear  types  of  representation.  More  Impor¬ 
tant,  the  simple  bilinear  model  Includes  no  accurate  account  of  the  strain  rate 
effect  on  the  material  properties.  We  Investigated  this  possibility  somewhat  by 
Increasing  the  yield  stress  until  finite  element  displacements  agreed  with  the 
experimentally  measured  values.  In  some  cases,  the  required  amount  of  yield  stress 
Increase  seemed  reasonable.  However,  the  bilinear  nature  of  the  stress-strain  curve 
and  lack  of  strain  rate  alteration  In  the  hardening  portion  remain  questionable 
areas.  The  effect  of  strain  rate  on  the  hardening  modulus  could  be  just  as  signifi¬ 
cant  as  yield  stress  In  predicting  displacement  when  plastic  zones  become  large.  If 
the  hardening  modulus  is  increased  only  slightly,  a  disproportionate  decrease  in  dis¬ 
placement  could  occur.  Investigation  of  this  phenomenon  Is  worthy  of  future  study, 
particularly  when  the  material  is  better  characterized. 

A  sizable  amount  of  current  research  is  concerned  with  development  of  accu¬ 
rate  finite  element  constitutive  relationships.  The  finite  element  code  ABAQUS 
offers  much  promise  in  this  area;  it  has  both  exponential  hardening  material  descrip¬ 
tion  and  Malvern  strain  rate  sensitivity  capabilities  built  in.^^  Had  this  code  been 
available  at  the  start  of  our  work,  it  might  have  been  used  for  the  analyses.  At 
this  writing,  consideration  Is  being  given  to  reanalyzing  some  of  the  specimens  in 
this  test  program  using  ABAQUS  and  a  presumably  better  material  description  that 
includes  rate  effects.  The  overall  lack  of  material  property  data  available  for 
HY-80  steel  tested  at  high  strain  rates  would  be  a  hindrance,  however.  This  lack  of 
data  makes  the  validity  of  the  derived  property  constants  that  are  required  question¬ 
able.  Therefore,  one  can  only  presume  that  having  this  strain  rate  capability  is 
truly  advantageous.  This  question  points  out  the  need  for  a  large  scale  material 
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characterization  program  for  the  naval  HY  steels  much  like  that  already  performed 
by  the  nuclear  industry  on  A533B  reactor  pressure  vessel  steel. 

A  third  possible  explanation  for  the  discrepancies  between  the  finite  element 
results  and  the  measured  test  results  is  Inaccuracy  in  the  finite  element  idealiza¬ 
tion.  As  mentioned  in  the  last  section,  inadequacies  in  the  side  groove  model  appear 
to  be  a  possible  cause  for  the  large  overpredlctlon  of  displacements  that  occurred 
in  the  analyses  of  the  more  severely  loaded  specimens,  and  particularly  for  the  speci¬ 
men  having  V-shaped  side  grooves.  Another  possible  Inaccuracy  in  the  finite  element 
model  is  the  manner  of  load  application  at  the  bolt  connection.  As  modeled  in  the 
finite  element  idealization  (Figure  22),  the  load  history  is  applied  as  a  concen¬ 
trated  load  at  one  point  in  the  specimen.  This  situation  would  correspond  to  a 
static  compact  specimen  test  with  undersized  load  pins.  For  this  dynamic  test  situa¬ 
tion,  tightly  fitting  bolts  were  used  to  transfer  the  load  in  order  to  remove  as 
much  of  the  "slack”  in  the  system  as  possible.  Considerable  effort  was  made  to 
allow  the  specimen  to  rotate  freely  about  these  bolts.  The  bolts  were  lubricated 
where  they  Interfaced  with  the  specimen  and  were  tightened  to  the  point  where  they 
could  still  spin  freely.  Despite  these  precautions,  there  may  have  been  some  applied 
moment  at  the  load  holes,  and  the  point  load  used  in  the  finite  element  analysis  may 
not  have  been  a  good  representation  of  the  distributed  loading  applied  by  the  snugly 
fitting  bolts.  This  is  another  area  for  future  investigation.  Idealization  changes 
in  the  load  pin  region  and  the  use  of  distributed  or  multiple  point  loads  may  cause 
unexpectedly  large  changes  in  the  displacement  response  of  the  specimen. 

The  finite  element  analyses  were  extremely  sensitive  to  crack  length.  This  was 
among  the  most  enlightening  results  of  the  various  finite  element  analyses  and  was 
in  fact  stumbled  upon  quite  accidently.  During  the  analysis  of  specimen  JC-27,  the 
crack  length  was  input  as  1.51  in.  Instead  of  the  actual  measured  value  of  1.50  in. 


When  this  was  corrected,  and  the  analysis  rerun  with  a  crack  depth  of  1.50  in.,  the 
result  was  a  40%  smaller  peak  load  point  displacement.  The  load  point  displacement 
histories  for  the  two  analyses  are  shown  along  with  the  measured  experimental  value 
in  Figure  36.  Both  analyses  were  conducted  using  a  yield  stress  20%  above  the  static 
value  and  the  same  elastic  and  hardening  moduli.  The  extent  of  the  change  in  dis¬ 
placement  response  was  quite  surprising.  To  further  investigate  this  effect,  another 
specimen  analysis  was  performed. 

Specimen  JC-23  was  more  deeply  cracked,  and  had  been  less  severely  loaded,  than 
specimen  JC-27.  We  chose  to  analyze  it  with  different  crack  depths  to  see  if  a 
change  in  the  load  point  displacement  as  large  as  that  found  in  specimen  JC-27  would 
occur  in  this  case.  Shown  in  Figure  37  are  the  LPD  histories  obtained  for  crack 
depths  of  1.58  in.  and  1.59  in.  As  in  specimen  JC-27,  this  small  change  in  crack 
depth  resulted  in  roughly  a  40%  difference  in  peak  displacement.  This  similar  result 
for  a  substantially  different  test  seems  to  Indicate  that  the  deeply  cracked  nature 
of  the  specimens  in  this  test  program  (a/W  from  0.725  to  0.79)  is  the  major  contrib¬ 
uting  factor  to  the  crack  depth  sensitivity. 

A  question  remains  as  to  whether  this  crack  depth  sensitivity  occurs  in  the 
actual  specimen  to  the  degree  in  which  it  occurs  in  the  finite  element  analyses. 

The  elastic  compliance's  versus  the  crack  depth  for  the  side  grooved  specimens  is 
plotted  in  Figure  38.  As  shown  in  this  figure,  the  a/W  values  for  the  specimens 
used  in  this  program  fall  in  a  very  rapidly  changing  section  of  the  curve.  Despite 
this,  even  in  the  most  deeply  precracked  specimens,  the  elastic  displacement  response 
would  only  vary  about  5%  for  a  crack  depth  change  of  0.01  in.  The  rest  of  the  change 
in  displacement  must  be  due  to  the  large  amount  of  plasticity  occurring  in  the  rela¬ 
tively  small  remaining  ligament  section.  As  shown  in  the  two  test  cases,  the  rela¬ 
tive  effect  due  to  change  in  crack  depth  appears  to  be  essentially  independent  of 
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yield  stress.  The  effect  of  variation  in  hardening  in  the  material  models  has  not 
been  Investigated  and  remains  a  questionable  area. 

The  ramifications  of  this  crack  depth  sensitivity  could  be  significant.  One 
must  question  the  finite  element  solution  for  any  very  deeply  cracked  specimen  for 
which  the  measured  crack  depth  Is  averaged  to  the  nearest  0.01  In.  and  input  In  the 
analysis.  These  results  Indicate  that  more  accurate  crack  depth  measurement  Is 
needed.  The  effect  of  curved  crack  fronts  also  appears  to  be  a  potential  source  of 
large  error  when  Idealized  as  a  planar  crack  front  In  a  two-dimensional  analysis. 
This  crack  sensitivity  could  have  been  a  contributing  factor  to  the  discrepancies 
between  the  experimental  and  finite  element  results  presented  earlier.  A  slight 
difference  between  the  actual  and  idealized  crack  depths  could  fully  explain  the 
overly  flexible  finite  element  solutions  obtained  using  what  are  considered  more 
reasonable  yield  stress  values.  Obviously,  this  topic  deserves  more  thorough  in¬ 
vestigation.  The  effect  of  static  versus  dynamic  loading  on  this  crack  depth  sensi¬ 
tivity  was  not  investigated.  More  complete  analyses  of  material  property  effects 
are  also  required.  A  great  deal  of  finite  element  work  has  been  performed  on  the 
analysis  of  compact  tension  specimens  (primarily  under  static  load)  and  good  agree¬ 
ment  between  experiment  and  analysis  has  often  been  obtained.  Why  then,  is  this 
extreme  crack  tip  sensitivity  observed  in  the  dynamic  analyses  performed  for  this 
work?  Possibly,  it  results  from  the  various  circumstances  involved  such  as  dynamic 
load,  deeply  precracked  specimens,  type  of  analysis,  finite  element  Idealization, 
and  material  properties.  Or,  possibly,  this  phenomenon  could  be  a  less  obvious 
source  of  error  in  other  work  involving  less  deeply  cracked  and  statically  loaded 
specimens  having  large  scale  plasticity. 
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SUMMARY  AND  CONCLUSIONS 


A  coitibined  experimental  and  analytical  research  program  has  been  described  that 
was  aimed  at  extending  a  static  nonlinear  fracture  criterion,  the  J-lntegral,  to 
inelastic  dynamic  fracture.  Within  this  program,  an  experimental  technique  was 
developed  by  which  standard  compact  tension  specimens  could  be  tested  dynam¬ 
ically,  and  Instrumentation  techniques  were  developed  which  permitted  measurement  of 
fracture  events  that  occurred  within  only  milliseconds.  Concurrently,  nonlinear, 
dynamic  finite  element  analyses  were  used  in  guiding  the  experimental  work  and  In 
interpreting  the  experimental  results.  The  work  described  was  a  pioneering  effort 
spanning  4  years.  We  encountered  many  problems,  particularly  during  the  experimental 
phase,  and  our  methods  of  dealing  with  these  problems  have  been  documented  in  detail. 

The  material  chosen  for  the  experimental  program  was  1 -in. -thick  HY-80  baseplate 
because  of  Its  known  toughness  and  widespread  use.  Machined,  precracked  compact 
tension  specimens  were  dynamically  loaded  by  a  weight  dropped  onto  a  special  test 
fixture  designed  to  load  In  the  opening  mode.  A  total  of  26  steel  specimens  were 
dynamically  tested  (some  repeatedly)  to  develop  and  prove  required  instrumentation 
and  to  measure  an  apparent  critical  value  of  the  J-integral.  Time  to  fracture  was 
of  the  order  of  1  ms  or  more.  Three  additional  specimens  were  tested  under  standard 
conditions  to  establish  the  static  value  of  Jjc* 

The  experimental  program  required  that  load  and  load  point  displacement  (LPD) 
be  measured  on  the  specimens  as  a  function  of  time.  Due  to  the  rapidity  of  the 
loading  and  deformation,  load  point  displacement  was  particularly  difficult  to  mea¬ 
sure.  The  best  system  devised  consisted  of  a  t r; nsmi t t Ing-r ece I vir g  photo-optic 
displacement  sensor  specially  calibrated  to  accoi nt  for  rotation  of  the  specimen 
about  Its  apparent  hinge  point.  This  gave  an  I.PD  accuracy  to  within  10%.  Load  was 
more  easily  measured  using  conventional  strain  gages  mounted  on  the  loading  clevis. 


To  determine  the  critical  value  of  the  J-lntegral,  It  was  also  necessary  to 
observe  the  time,  load,  or  LPD  at  which  fracture  Initiation  occurred.  Various  tech¬ 
niques  were  tried,  but  even  the  most  promising  (DC  potential  drop)  could  not  be 
calibrated  to  reliably  Indicate  fracture  Initiation.  Because  of  the  large  plastic 
strains  that  develop  near  the  crack  tip  before  fracture,  strain  gages  and  "Krak" 
gages  were  also  Ineffective;  they  came  unbonded  from  the  specimen  surface  before  the 
start  of  crack  growth.  As  a  result,  a  multiple  specimen  approach  was  adopted  in 
which  specimens  were  Incrementally  loaded  up  to  and  just  beyond  the  point  of  fracture 
Initiation.  This  yielded  an  "apparent  Jic"  value,  as  calculated  by  conventional 
static  formulas,  of  about  1000  In.-lb/ln.^  for  the  HY-80  material  under  test. 

This  compares  to  an  average  static  value  of  1110  ln.-lb/ln.2. 

The  term  "apparent  Jic"  is  used  above  for  good  reason.  Plots  of  load  versus 
LPD  will  differ,  at  a  given  value  of  load  or  displacement,  under  static  and  dynamic 
conditions.  The  additional  load  on  a  dynamically  tested  specimen  at  the  same  value 
of  static  displacement  will  depend  on  the  rate  of  loading,  and  this  will  influence 
the  value  of  the  J-integral  calculated  from  the  experimental  record.  Nonetheless, 
direct  use  of  the  dynamic  records  has  been  adopted  in  the  past^^>3^  and  is  adopted 
here.  As  a  consequence,  the  exact  meaning  of  the  fracture  parameters  calculated 
herein  is  presently  unknown.  This  question  is  currently  under  active  investigation 
by  a  number  of  well-known  researchers. 3^. 35  Despite  the  uncertainty,  however,  the 
comparative  fracture  performance  of  different  materials  can  be  evaluated  with  the 
present  experimental  setup. 

The  finite  element  computations  were  an  important  part  of  this  work,  and  led  us 
to  correctly  question  the  accuracy  of  some  of  the  early  experimental  data.  Another 
motive  for  undertaking  these  calculations  resulted  from  the  experimental  difficulty 
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of  detecting  crack  Initiation.  The  thinking  was  that  If  one  can  closely  model  the 
experiment  numerically,  then  any  deviation  between  the  experiment  and  the  numerical 
model  Is  due  to  crack  extension,  thus  pinning  down  that  event. 

A  complicating  factor  In  the  two-dimensional  analysis  of  these  specimens  was 
the  presence  of  side  grooves  designed  to  suppress  the  formation  of  shear  lips.  These 
were  modeled  by  reducing  the  thickness  of  the  elements  In  the  side  groove  region  and 
by  reducing  Young's  modulus  In  this  region  In  order  to  bring  the  finite  element 
static-elastic  compliance  Into  agreement  with  accepted  calibration  formulas.  The 
suitability  of  this  model  In  the  dynamic,  elastic-plastic  range,  however,  can  be 
established  only  by  comparison  with  fully  three-dimensional  solutions.  The  expense 
of  such  solutions  precluded  this  check. 

Another  complicating  question  Is  whether  a  model  of  plane  stress  or  plane  strain 
is  most  appropriate  for  what  Is,  In  reality,  a  fully  three-dimensional  problem.  Near 
the  crack  tip,  plane  strain  Is  probably  more  appropriate,  with  plane  stress  preferred 
elsewhere.  Neither  assumption  Is  exactly  correct,  but  a  plane  stress  model  was 
assumed  for  all  computations  in  this  study. 

Finally,  It  Is  very  Important  to  note  that  the  accuracy  of  the  finite  element 
models  was  quite  sensitive  to  assumptions  for  the  stress-strain  curve  and  the  Initial 
crack  depth.  We  found  that  a  5%  Increase  In  assumed  yield  stress  could  decrease  the 
maximum  load  point  displacement  by  about  40%.  This  indicates  that  the  effect  of 
strain  rate  on  yield  stress  should  be  Incorporated  into  the  analysis.  We  also  found 
that,  for  these  deeply  cracked  specimens,  an  error  in  the  crack  depth  of  less  than 
1%  (i.e.,  an  error  of  O.Ol  In.)  also  could  change  the  peak  load  point  displacement 
by  40%.  These  unwanted  sensitivities  arise  because  of  the  extensive  plasticity  that 
develops  In  the  specir  2ns  and  because  of  the  rapid  change  in  (nonlinear)  specimen 
compliance  for  small  changes  in  crack  depth.  Because  the  crack  front  Is  rarely 
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perfectly  straight,  measuring  crack  depth  precisely  Is  difficult,  which  makes  the 
sensitivity  of  peak  displacement  to  this  parameter  particularly  disturbing.  This 
result  Implies  that  deeply  cracked,  highly  nonlinear  compact  tension  specimens  are 
not  the  best  choice  for  characterizing  Inelastic  dynamic  fracture  In  very  ductile 
metals. 

The  purpose  of  this  work  was  not  to  characterize  materials,  however,  although 
the  experimental  techniques  and  analysis  methods  are  certain  to  be  of  future  use. 
Instead,  we  sought  to  determine  If  the  fracture  response  of  HY  steels  Is  signifi¬ 
cantly  different  under  static  and  dynamic  conditions.  If  so,  the  difference  would 
have  to  be  accounted  for  In  the  Inelastic  fracture  analysis  of  full  scale  structures. 
This  work  gives  the  preliminary  Indication  that  at  nominal  strain  rates  of  Interest, 
the  difference  Is  not  great  enough  to  be  significant.  However,  this  conclusion  must 
be  qualified  by  noting  that  It  Is  based  on  the  assumption  that  the  J-llke  fracture 
criterion  derived  here  has  roughly  the  same  meaning  as  the  J-lntegral  under  static 
conditions. 
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FURTHER  INSTRUMENTATION  DETAILS 
MEASUREMENT  OF  LOAD  POINT  DISPLACEMENT 

The  sensor  Initially  selected  was  a  displacement  measuring  system  series  KD-2310 
manufactured  by  Kaman  Sciences.  Its  body  was  1  1/8-in.  long  1/4-ln.  diameter  and 
was  mounted  on  the  end  of  a  small  diameter  cable.  It  was  complete  with  a  control 
module  having  provisions  for  adjusting  the  gain  and  linearity.  According  to  specifi¬ 
cation,  it  had  a  frequency  response  from  0  to  50  kHz  and  a  linear  response  to  55  mil 
(1  mil  =  0.001  in.).  It  was  a  noncontact  type  transducer  that  operated  on  the  eddy- 
current  loss  principal  and  could  be  used  with  any  conducting  material.  Another 
advantage  was  that  it  Indicated  the  average  gap  over  the  face  of  the  probe.  Thus, 
providing  and  maintaining  a  true  parallel  target  surface  was  not  essential. 

One  disadvantage  of  this  sensor  was  that  metal  located  to  the  side  and  rear  of 
its  face  could  Influence  the  sensitivity.  To  reduce  this  effect,  the  sensor  type 
having  a  stainless  steel  jacket  was  selected  and  calibrated  in  a  mounting  configur¬ 
ation  similar  to  that  used  in  the  test,  as  recommended  by  the  manufacturer.  A 
micrometer  fixture  was  used  for  calibrations. 

The  sensor  was  threaded  through  the  top  half  of  the  specimen  and  targeted  on 
the  lower  half;  initially  a  spacer  was  bonded  in  the  notch  and  the  sensor  was  located 
flush  with  the  bottom  of  the  hole.  This  would  protect  the  sensor  If  the  notch  closed 
during  latter  phases  of  the  test.  However,  this  problem  did  not  develop,  so  the 
spacer  was  eliminated.  Then  the  sensor  was  referenced  directly  on  the  bottom  face 
of  the  notch.  This  alternate  procedure  helped  eliminate  the  Influence  of  surrounding 
metal,  which  we  believed  was  causing  sensitivity  shifts. 
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After  some  actual  tests,  certain  problems  with  the  initially  selected  sensor 
began  to  emerge.  The  major  problem  was  that,  at  fracture  initiation,  the  specimen 
appeared  to  open  more  than  anticipated,  exceeding  the  linear  travel  of  the  sensor. 

As  the  testing  progressed,  this  was  confirmed  by  comparison  with  other  instrumenta¬ 
tion,  and  also  by  physical  measurements  taken  after  the  test.  In  addition,  the 
system  seemed  sensitive  to  movements  of  the  sensor  cable.  This  Influence  was  reduced 
by  affixing  the  cable  and  module  to  a  long  board,  leaving  about  2  ft  free  to  extend 
to  the  sensor.  Another  disturbing  problem  was  the  shifting  of  the  sensitivity  for  no 
apparent  reason.  Figure  A.l  shows  the  setup  and  sensitivity  for  the  sensor  before 
using  it  in  a  test  and  the  calibration  Immediately  afterwards.  Correlation  was  good 
up  to  about  25  mil,  but  beyond  this,  the  sensitivity  after  the  test  increased 
nonllnearly.  In  comparison,  the  sensitivity  before  the  test  was  linear  up  to  approx¬ 
imately  55  mil.  This  particular  test  was  at  a  moderate  drop  and,  by  using  the 
pretest  sensitivity,  the  dlplacement  was  52  mil.  However,  using  the  latter  sensi¬ 
tivity,  it  would  be  43.5  mil.  This  represents  a  variation  of  roughly  10%  from  the 
mean.  This  discrepancy  was  more  extreme  than  was  usually  noticed,  but  it  certainly 
diminished  the  confidence  in  the  sensor.  The  combination  of  cable  shifting  and 
influence  of  surrounding  metal,  despite  the  precautions  taken,  probably  caused  this 
change.  This  same  type  sensor  was  used  at  position  D2  and  eventually  at  D2'  (see 
Figure  7).  The  accuracy  of  the  measurement  at  either  of  these  positions  was  not 
critical,  however,  as  will  be  discussed. 

Because  of  the  erratic  performance  and  limitations  of  the  eddy-current  sensor. 
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it  became  apparent  that  an  alternative  type  was  needed  to  measure  the  LPD.  Thus,  a 
basically  different  system  was  obtained  on  a  temporary  basis  for  evaluation,  and 
eventually  this  system  was  purchased.  As  described  in  the  text,  this  system  was  a 


Fotonic  sensor,  series  KD320  with  a  KDP109  CTI  probe,  manufactured  by  Mechanical 


Technology,  Inc.  The  sensing  probe  contains  two  sets  of  fiber  optic  filaments.  The 
system  had  a  frequency  response  from  DC  to  125  kHz  according  to  specifications.  The 
probe  was  flexible  except  for  a  6-ln.  steel  jacket  at  the  sensing  end.  The  sensor 
had  to  approach  from  the  side  and  enter  the  specimen  from  the  top.  For  the  first 
attempt  the  probe  was  custom  configured  to  Include  a  right-angle  bend  near  the  end 
so  that  with  some  care  It  could  be  adapted  to  the  specimen.  This  arrangement  was 
too  susceptible  to  damage  and,  after  consultation  with  the  manufacturer.  It  was 
determined  that  the  rigid  section  could  be  shortened  and,  after  some  machining  opera¬ 
tions  to  provide  more  clearance,  the  flexible  cable  could  take  the  necessary  bend. 

By  using  an  adaptor  this  sensor  was  mounted  in  the  same  hole;  it  was  targeted  on  the 
bottom  of  the  notch. 

During  a  test  the  specimen  can  be  evisioned  as  hinging  in  the  ligament  section 
with  the  upper  and  lower  halves  counterrotating  with  respect  to  each  other.  For  a 
realistic  test  the  included  angle  from  this  rotation  can  approach  4°.  This  means 
that  target  face  does  not  remain  perpendicular  to  the  sensor  and  the  impinging 
light  from  the  surface  would  tend  to  be  scattered  during  the  test  and  influence  the 
sensitivity.  This  factor,  of  course,  must  be  considered  during  calibration  and 
adjustment  of  the  measuring  system  prior  to  use.  The  micrometer  fixture  again  was 
used  to  adjust  and  calibrate  the  system,  although  a  modification  was  incorporated 
to  simulate  this  rotating  effect  (Figure  A. 2).  Note  the  hinge  Installed  in  this 
fixture,  which  pivots  with  motion  of  the  micrometer  head.  The  length  L  corresponds 
to  the  distance  between  the  sensor  and  the  estimated  hinge  point  in  the  specimen 
when  loaded  as  computed  from  Information  in  Reference  21.  Figure  A. 3  shows  a  typical 
calibration  curve.  Note  there  are  front  slope  and  back  slope  operating  regions. 
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The  rear  slope  would  be  used  in  situations  such  as  these  dynamic  tests,  where  larger 
displacements  are  Involved. 

Different  probe  diameters  and  arrangements  for  the  optical  fibers  Influence  the 
operating  ranges.  For  this  application  It  was  necessary  to  use  the  smallest  avail¬ 
able  probe  (0.109-in.  diameter)  and,  as  a  result,  the  available  displacement  range 
was  sacrificed.  The  probe  had  a  concentric  fiber  distribution  where  the  transmitting 
fiber  formed  the  inner  core  and  the  perceiving  fibers  the  outer  ring.  This  arrange¬ 
ment  was  chosen  because  It  was  less  sensitive  to  the  nonperpendicular  target  face. 
Figure  A. 4  shows  a  typical  sensitivity  curve  for  this  sensor  where  the  target  face 
remains  perpendicular  to  the  sensor  and  a  second  curve  where  the  target  face  Is 
initially  perpendicular  and  then  rotates  to  match  the  situation  In  the  test  specimen. 
Note  that  the  rotating  target  face  affects  the  relationship,  especially  the  sensitiv¬ 
ity  in  the  back  slope  region.  Also,  the  sensitivity  peaks  at  a  different  displacement 
for  the  two  target  face  arrangements.  To  operate  in  the  linear  range  of  the  back 
face  slope  region,  an  initial  gap  of  70  to  90  mil  is  required.  Because  of  this  and 
also  because  In  the  specimen  the  target  face  is  Initially  perpendicular  at  the  time 
of  loading,  an  additional  step  is  required  to  obtain  a  valid  sensitivity.  Repeated 
calibrations  for  different  initial  gaps  indicated  that  80'  mil  gave  optimum  useful 
range  (Figure  A. 5).  The  sensitivity  is  calculated,  and  within  the  *5%  bounds,  this 
calculated  figure  is  valid  for  more  than  about  140  mil  displacement.  In  practice, 
the  sensor  was  therefore  installed  with  an  80-mil  gap  and,  through  a  control  on  the 
instrument,  the  light  intensity  was  adjusted  to  give  the  corresponding  output  volt¬ 
age.  This  adjustment  was  necessary  to  compensate  for  different  remlssivlty  between 
target  surfaces.  To  minimize  this  effect,  a  highly  reflective  shim  material  was 
always  bonded  to  the  target  face  both  during  the  calibration  and  the  testing 
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To  evaluate  the  Fotonic  sensor  under  actual  tests,  a  third  tranducer  was  added 
to  the  specimen  and  was  also  recorded  during  several  runs.  This  additional  sensor 
was  a  linear  variable  differential  transformer,  LVDT;  it  was  epoxled  to  the  front 
face  and  also  measured  the  displacement.  To  compare  results,  we  again  assumed  that 
the  specimen  hinges  in  the  ligament  about  a  calculated  center  of  rotation.  Reference 
21  discusses  how  this  rotation  point  can  be  determined.  It  states  that  for  an 
elastic  strain  at  the  crack  tip  this  rotation  would  occur  at  a  point  10%  through  the 
ligament  from  the  crack  tip  and  for  a  large  amount  of  plasticity  at  roughly  33%  of 
the  ligament  length.  The  displacement  magnitudes  Indicated  by  each  sensor  would  be 
proportioned  to  its  distance  from  the  center  of  rotation  along  a  base  line.  In 
Figure  A. 6  the  results  from  these  three  sensors  are  plotted  at  corresponding  loca¬ 
tions  for  two  different  static  loads  on  one  specimen  and  for  an  actual  drop  test  on 
a  second  specimen.  Note  that  for  the  static  test  the  correlation  was  very  good,  but 
for  the  impact  test  the  Fotonic  sensor  was  8%  higher  than  expected.  In  Figure  A. 6b 
the  displacements  observed  after  the  test  (set  values)  are  shown  and  the  correlation 
between  the  three  sensors  is  good.  This  certifies  that  the  sensitivities  for  the 
sensors  were  correct.  Thus  this  8%  disparity  in  the  test  values  is  difficult  to 
explain.  It  cannot  be  attributed  to  nonlinearity  because  the  calibration  curves 
show  a  constant  relationship  over  this  range.  Based  on  this  study,  as  well  as  many 
similar  cross  checks  made  between  the  displacements  thoughout  the  test,  it  appears 
that  the  LPD  measurement  using  the  Fotonic  sensor  can  be  considered  accurate  to 
within  10%. 

CRACK  INITIATION 

The  next  section  is  devoted  to  the  various  procedures  attempted  to  detect  crack 
initiation.  A  crack  would  be  expected  to  initiate  at  mldthlckness  and  this  was  a 


DISPLACEMENT  (mil)  DISPLACEMENT  (mil) 


LOCATIONS  ON  SPECIMEN 


DISPLACEM^T  SENSOR  TYPE 

LVDT-LINEAR  VARIABLE  DIFFERENTIAL  TRANSFORMER 
D1-FOTONIC  SENSOR 

© 

D2  -EDDY  CURRENT  SENSOR 

f 

LVDT 

toi  j 
© 

Figure  A. 6a  -  Displacement  Measured  Under  Static  Load 
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consideration  in  the  original  plans.  At  first  we  used  an  eddy-current  sensor  iden¬ 
tical  to  the  one  discussed  earlier  in  connection  with  the  LPD  measurement.  It  was 
threaded  into  a  blind  hole  from  the  top  of  the  specimen  Into  the  ligament  section. 
This  sensor,  designated  D2 ,  targeted  the  bottom  of  the  hole  ahead  of  the  crack  front 
Figure  A. 7  shows  this  location  and  also  the  method  by  which  we  hoped  to  detect  the 
crack  initiation.  This  ligament  displacement,  LD,  and  the  LPD  are  taken  from  the 
record  and  plotted  at  corresponding  time  Increments.  A  sudden  change  In  the  slope 
would  indicate  the  onset  of  fracture.  After  several  tests  we  concluded  that  the 
ligament  gage  was  insensitive  to  crack  Initiation  and  any  significant  response  oc¬ 
curred  only  when  the  crack  actually  reached  the  sensor;  this  was  evident  from  inspec 
tlon  of  the  specimen,  from  time  relations,  as  well  as  from  the  magnitude  of  the 
measurement.  Figure  A. 8  is  a  continuous  record  of  the  applied  load,  the  load  point 
displacement,  and  the  ligament  displacement  recorded  during  a  test.  Note  that 
there  is  no  signal  on  the  LD  trace  until  nearly  5  tns,  after  which  the  displacement 
steadily  increases  and  peaks  at  roughly  25  rail  in  another  5  ms.  Visual  inspection 
revealed  that  the  crack  in  the  specimen  obviously  had  extended  beyond  the  sensor, 
accounting  for  this  large  displacement.  Attempts  to  shift  the  ligament  gage  closer 
to  the  crack  tip  for  more  sensitivity  were  unsuccessful.  There  was  concern  that  the 
sensor  mounting  hole,  if  located  closer  to  the  crack  tip,  would  cause  premature 
fracture,  and  therefore  the  approach  was  changed.  Instead  of  being  placed  in  the 
ligament  the  sensor  was  located  at  the  root  of  the  machined  notch  and  targeted  on 
the  Lower  face.  Figure  A. 9  shows  results  for  this  configuration  from  two  tests 
where  it  was  evident  from  examination  that  both  specimens  had  cracked.  Slope  did 
not  change  definitively  even  though  fracture  occurred. 
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Figure  A. 8  -  Attempt  to  Measure  Crack  Initiation 
with  Sensor  Located  in  Ligament 
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Figure  A. 9  -  Attempts  to  Measure  Crack  Initiation 
with  Sensor  Treated  in  Notch 
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At  this  point  the  emphasis  was  shifted  to  evaluation  of  a  number  of  different 
procedures;  they  Included  using  crack  detection  and  crack  propagation  gages, 
moire  fringes,  and  finally  the  potential  difference  (PD)  method. 

A  number  of  different  types  of  gages  were  bonded, to  the  sides  of  the  specimen 
In  the  vicinity  of  the  crack,  as  shown  In  Figure  10  of  the  main  text.  At  the  very 
beginning  of  this  test  program,  a  bonded  gage  was  tried  briefly  with  the  Idea  of 
providing  supplemental  data;  see  Figure  10a.  This  was  a  crack  propagation  gage  and 
consisted  of  a  pattern  of  resistance  strands  connected  In  parallel.  It  was  bonded 
with  the  strands  running  normal  to  the  crack  front.  Progression  of  the  crack 
causes  successive  Interruptions  of  the  strands,  accompanied  by  an  Increase  In  resist¬ 
ance.  With  proper  circuitry,  this  change  Is  recorded  as  discrete  steps  from  which 
the  crack  front  can  be  tracked.  After  tests  on  the  first  specimen  It  was  observed 
that  strands  had  broken  ahead  of  where  the  crack  terminated  and  there  was  no  further 
use  of  this  gage  type. 

When  work  was  resumed  using  bondable  gages,  several  specimens  were  prepared 
with  a  small  strain  gage  Installed  on  the  side  faces  In  the  vicinity  of  the  crack 
tip.  This  was  attempted  both  with  the  side  surfaces  flat  and  with  the  wide  side 
grooves.  The  thought  was  that  with  this  gage  In  close  proximity  to  the  crack  tip, 
there  would  be  a  strain  release  at  crack  Initiation  which  would  be  evident  on  the 
record.  Figures  10b  and  10c  show  the  layout  for  these  bonded  strain  gages.  Another 
procedure  using  "Krak"  gages  was  also  tried  In  an  effort  to  detect  crack  Initiation. 
Figure  lOd  shows  tha  Installation  arrangement  where  the  gage  was  placed  directly 
over  the  crack  tip  to  follow  the  fracture  propagation  In  the  specimen.  This  gage 
consisted  of  conducting  foil  on  a  thin  backing  material.  It  was  bonded  using 
conventional  material  and  techniques  and  It  operated  In  a  Wheatstone  bridge  fashion 
using  a  constant  current  excitation. 


It  soon  became  evident  that  bonded  gages  were  not  satisfactory;  they  almost 
Invariably  came  loose  during  the  test.  The  data  showed  this  problem,  and  some  gages 
were  actually  visibly  separated  from  the  specimen.  Figure  A. 10  shows  records  from  a 
strain  gage  and  a  "Krak“  gage  located  on  opposite  sides  of  the  specimen  along  with 
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the  applied  load  taken  during  a  test.  In  this  parlcular  test  the  specimen  did  not 
fracture,  although  there  Is  a  sudden  jump  In  each  trace.  Note  also  that  for  the 
crack  gage  the  signal  occurs  early  In  the  event  when  the  load  Is  roughly  one-half 
the  peak  value,  nearly  1  ms  ahead  of  the  strain  gage.  Both  gages  obviously  came 
unbonded.  In  some  cases  there  was  evidence  that  the  strain  gages  did  not  always 
come  loose.  For  example.  Figure  A. 11  shows  results  for  a  strain  gage  bonded  In  each 
groove.  This  specimen  did  not  crack  and  there  Is  no  sudden  shift  In  the  signal 
(perfect  correlation).  However,  Figure  A. 12  Is  a  similarly  prepared  specimen  where 
crack  Initiation  occurred,  and  this  shows  an  entirely  different  problem.  One  gage 
showed  a  sudden  unloading  while  the  other  did  not  and  was  obviously  Insensitive  to 
crack  Initiation.  Thus  It  was  apparent  from  this  work  that  bonded  strain  gages 
were  not  a  solution.  They  were  either  Insensitive  or  they  came  unbonded  and  gave 
meaningless  data. 

A  number  of  other  problems  were  noticed  with  the  "Krak”  gages,  several  of  which 
also  substantiated  this  bonding  difficulty.  These  gages  were  Installed  after  the 
specimen  had  been  precracked,  and  therefore  the  specimen  had  to  be  exercised  for 
several  cycles  to  produce  a  corresponding  crack  In  the  gage.  In  a  number  of  cases 
It  was  apparent  that  the  crack  ^n  the  gage  was  not  following  the  crack  In  the 
specimen,  and  In  fact  the  gage  had  come  loose  near  the  crack  tip.  In  one  case  the 
gage  came  unbonded  during  the  Initial  loading  without  a  crack  even  starting.  Vfhen 
bonding  problems  were  detected  the  gages  were  replaced  until  Installation  seemed  to 
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Figure  A. 11  -  Example  Showing  Good  Correlation  Between 
Crack  Detection  Gages 
(No  Crack  Extension) 
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Figure  A. 12  -  Example  Showing  Discrepency  Between  Crack  Detection 
Gages  (Crack  Extension  Measured  in  Specimen) 
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be  satisfactory.  For  Installation  the  gages  were  fitted  Into  the  groove,  which 
seemed  to  compound  the  problem,  as  there  was  tendency  for  the  gage  to  separate  over 
a  wide  area.  Another  problem  was  the  tendency  during  Installation  for  the  gage  to 
open  and  lose  continuity  In  the  narrow  sections  where  the  gage  wraps  over  the  outside 
corner  of  the  specimen,  even  though  rather  large  radii  had  been  ground  at  these 
locations.  A  smaller  gage  which  essentially  fit  into  the  groove  was  subsequently 
used  to  solve  this  problem.  The  original  gages  had  a  10  mm  useful  length.  The 
smaller  gage  had  a  5-mra  length,  which  was  adequate,  but  It  had  to  be  more  carefully 
placed.  Another  difficulty  was  that  once  the  crack  developed  and  the  load  was  re¬ 
moved,  the  gage  actually  closed  so  that  electical  conduction  was  restored  and  Indi¬ 
cated  no  crack  growth.  Figure  A.  13  shows  results  from  two  Instrumented  specimens 
that  had  previously  been  exercised  to  produce  cracks  in  the  gages.  The  gages  were 
monitored  as  the  load  was  applied;  note  that  a  load  of  roughly  800  lb  was  required 
before  the  gages  were  fully  opened  to  indicate  the  true  crack  length.  In  an  attempt 
to  offset  this  effect,  two  specimens  were  held  under  a  compressive  preload  while 
the  gages  were  Installed.  Based  on  compliance  relations,  this  preload  was  roughly 
1000  lb.  Still,  after  the  gages  had  been  exercised,  they  also  effectively  closed 
when  the  load  was  dropped.  A  high  gain  in  the  recording  channel  was  desired  in 
order  to  detect  crack  initiation,  but  this  was  not  possible  without  special  equipment 
to  suppress  this  significantly  larger  precrack  signal.  Therefore,  the  choice  was 
made  to  separate  the  crack  faces  by  raking  the  gage  with  a  small  scribe.  This  was 
done  with  the  aid  of  a  microscope  and  the  gage  was  monitored  so  that  the  final 

indication  observed  during  the  exercising  process  was  obtained.  This  procedure  was 
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tedious;  furthermore,  the  possibility  of  precisely  matching  the  crack  was  question¬ 
able.  In  the  end  it  was  apparent  that  crack  gages  were  not  a  viable  solution  for 
detecting  crack  initiation  in  this  application. 
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Figure  A. 13  -  Load  Required  to  Open  "Krak" 

Gage  After  Precracking 

The  moire  fringe  technique  was  used  briefly  with  the  idea  that  strains  as 
well  as  crack  initiation  could  be  detected.  This  method  uses  an  optical  effect 
caused  by  a  grid  that  is  distorted  along  with  the  surface  and  interacts  with  a  second 
stationary  grid.  In  this  application  one  grid  is  cemented  over  the  crack  tip  section 
of  the  specimen  and  the  second  grid  overlaid  and  supported  separately;  they  are 
referred  to  as  the  working  and  the  master  grid,  repsectively.  These  grids  consist 
of  unidirectional  lines;  as  the  specimen  distorts,  the  superposition  causes  inter¬ 
ference  and  moire  fringes  appear.  If  no  other  Influence  is  present,  these 


fringes  represent  a  map  of  the  strains  In  the  direction  perpendicular  to  the  lines 
In  the  master  grid.  The  magnitude  for  these  strains  can  be  determined  from 

e  =  p/d 

where  e  Is  the  strain,  P  Is  the  spacing  for  the  master  grid  and  d  the  distance  be¬ 
tween  observed  fringes.  The  spacing  for  the  master  grid  was  200  llnes/ln.  and  5Z 
was  taken  as  a  resonable  strain  at  the  crack  tip.  This  strain  would  produce  fringes 
0.01  In.  apart.  Thus  this  approach  seemed  suitable  for  measuring  the  plastic  strain 
around  a  crack  tip.  High  speed  16-mm  films  were  taken  of  these  fringes  during  some 
tests.  High  Intensity  flood  lamps  were  used,  and  by  forcing  the  developing  process, 
films  of  reasonable  clarity  were  obtained,  but  resolution  was  not  good  enough  to 
give  all  the  desired  Information.  Fringes  were  observed  but  they  were  related  to 
another  phenomenon  more  dominate  than  that  associated  with  tip  strains.  The  speci¬ 
men  halfs  tend  to  counterrotate  as  two  sections  due  to  the  bending  moment  at  the 
time  of  loading.  Fringe  patterns  consisting  of  a  series  of  contlnous  parallel  lines 
are  developed  as  a  result  of  this  rotation.  In  Figure  A. 14  are  two  enlargements 
obtained  from  the  16-mm  film  showing  the  zero  time  and  this  fully  developed  condi¬ 
tion.  The  distance  between  the  fringes,  d,  is  related  to  this  angle  of  rotation,  6, 
by 
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d  =  2(sln  0/2) 

A  typical  value  for  the  Included  angle  of  rotation  is  4°  and  where  P  =>  0.005  (200 
lines/ln.)  the  value  of  d  would  be  0.14  In.  This  Is  an  order  of  magnitude  greater 
than  the  separation  expected  at  the  crack  tip,  thus  overwhelming  the  fringe  pattern. 
A  time  study  was  made  of  this  angle  for  some  specimens.  Through  geometrical  re¬ 
lations,  the  angle  was  transformed  Into  displacement  for  correlation  with  the 


Figure  A. 14  -  Moire  Fringes  Produced  During  a  Severe  Specimen  Test 


measured  load  point  value;  see  Figure  A. 15  for  typical  results.  In  this  figure  the 


LPD  measured  by  the  instrumentation  is  compared  with  results  from  the  high  speed 
film  that  included  the  moire  grid  conversion  and  the  direct  measurements.  The 
values  from  the  instrumentation  and  directly  from  the  film  are  comparable,  consider¬ 
ing  that  the  magnitude  is  in  the  range  where  the  accuracy  of  the  eddy  current  sensor 
is  suspect  (see  Figure  A.l).  The  displacement  from  the  moire  fringe  pattern  is 
roughly  17%  low.  This  discrepancy  is  probably  due  to  the  Inaccuracy  of  reading  the 
spacing  of  the  fringe  because  of  poor  contrast  and  broad  lines  on  the  image  (see 
Figure  A. 14). 


Figure  A.  15  -  Comparison  of  Measurements  from  High  Speed  Film  and  LPD 
Sensor  for  Specimen  JC-15 


The  potential  difference,  PD,  was  the  final  method  tried  in  the  effort  to  detect 
crack  initiation.  This  technique  relies  upon  the  increase  of  resistance  as  the 
cross-sectional  area  decreases  as  in  the  case  for  crack  growth  in  the  CT  specimen. 

A  direct  current  is  passed  through  the  specimen  to  produce  a  potential  related  to 


this  uncracked  area.  Crack  initiation  would  of  course  be  indicated  when  the  poten¬ 
tial  exceeds  the  final  level  observed  during  precracking. 

In  Reference  22,  an  analytical  solution  was  derived  to  express  the  PD  through  a 
thin  bar  type  specimen  containing  a  small  crack  (a  machined  slot  in  this  case).  By 
solving  Laplace  equations  for  initial  conditions  and  using  conformal  mapping  proce¬ 
dures  the  potential  distribution  for  this  bar  is 
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where  V  is  the  potential  difference  between  two  points  on  the  specimen  center  line 
spaced  a  distance  Y  on  either  side  of  the  crack,  Sq  is  the  initial  crack  length 
and  a,  the  length  after  growth.  Solutions  have  been  expanded  to  cover  more  complex 
sections  and  results  have  been  confirmed  experimentally.  The  CT  specimen  is  too 
geometrically  complex  for  a  similar  approach  but  numerical  solutions  using  finite 
element  methods  have  been  successful  and  results  have  been  confirmed  experimentally 
(Reference  24).  In  the  experimental  work,  the  length,  a,  was  simulated  by  saw- 
cutting  the  specimen  in  incremental  steps  and  measuring  the  corresponding  potential 
change.  Reference  22  showed  that  the  attachment  position  for  the  probes  has  a  large 
effect  on  the  sensitivity  of  the  signal.  Following  these  guides,  the  potential 
probes  were  located  directly  across  the  plane  of  the  crack,  because  this  location 
produces  a  relatively  high  potential  change  with  crack  growth.  To  determine  the 
Influence  of  the  current  probes  on  sensitivity,  some  preliminary  work  was  done  where 
two  different  locations  on  the  top  and  bottom  surfaces  were  investigated  (see  Figure 
11).  Ring  lugs  were  soldered  to  the  ends  of  the  cables  and  connections  were  made 


with  screws  that  were  threaded  into  the  specimen.  A  block  diagram  of  the  circuit 
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used  In  this  testing  Is  shovm  In  Figure  A. 16.  The  PD  characteristic  was  observed 
during  the  fatigue  cracking  of  a  number  of  specimens.  Compliances  were  measured 
during  this  operation,  and  by  use  of  the  empirical  relations  described  in  Reference 
25  the  crack  growth  could  be  determined.  In  Figure  A. 17  the  voltage  observed  period¬ 
ically  at  peak  load  and  the  normalized  crack  growth  are  plotted  for  two  specimens, 
with  the  current  probes  located  at  the  two  different  positions  shown  in  Figure  11. 
Note  that  the  voltage  differed  by  a  factor  ranging  from  2.2  to  2.0  when  the  probes 
were  located  1-1/2  In.  rather  than  1/4  in.  from  the  back  edge  of  the  specimen.  On 
another  specimen,  which  was  already  precracked,  this  influence  was  determined  by 
shifting  these  connections  between  the  two  locations.  The  specimen  was  installed  in 
the  fixture  but  there  was  no  applied  load  at  the  time  of  the  readings.  The  voltage 
was  greater  by  a  factor  of  2.1  to  1,  which  agrees  with  this  earlier  observation.  In 
Figure  A. 18a  the  potential  difference  is  plotted  against  the  crack  growth  obtained 
during  precracking  where  the  current  probes  were  located  1-1/2  in.  from  the  back 
edge  of  the  specimen.  Both  parameters  have  been  normalized  to  the  initial  potential 
and  the  width  of  the  specimen,  respectively.  For  comparison,  similar  experimental 
results  taken  from  Reference  24  are  also  shown.  At  the  bottom  of  the  figure  the 
dimensions  for  these  two  specimens  are  compared.  Note  they  are  very  similar  except 
for  the  notch  depth,  and  this  accounts  for  the  two  different  a/W  realms  in  the  plot. 
Both  studies  show  similarity  In  chat  the  potential  varies  exponentially.  Where  the 
a/W  ratios  coincide,  you  might  expect  the  slopes  to  agree;  however,  note  that  the 
curve  from  the  present  data  is  somewhat  steeper.  This  may  be  attributed  to  a  number 
of  factors,  primarily  the  connection  points  for  the  Input  current  (compare  Figures 
11  and  A. 18b  for  details).  The  test  fixture  characteristics  can  also  affect  this 
sensitivity,  as  it  provides  a  parallel  electrical  path. 
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Figure  A. 18a  -  Comparison  of  Present  and  Referenced  Data 
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Figure  A. 18b  -  Comparison  of  Dimensions 


Figure  A. 18  -  Specimen  Potential  Drop  Characteristics  for  Present 

Test  and  Referenced  Results 


Another  was  that,  once  the  specimen  was  precracked  and  the  load  dropped,  the 
crack  faces  in  general  did  not  make  perfect  electrical  contact,  and  in  some  cases, 
there  was  little  difference  in  the  potential  between  the  peak  cyclic  load  and  no 
load.  Of  nine  specimens  tested,  eight  showed  a  drop  of  less  than  7%.  Halllday  and 
Beevers23  discussed  this  crack  face  contact  problem  and  inconsistencies  in  behavior 
during  cyclic  loads.  These  same  authors  caution  that  this  method  may  be  less  accu¬ 
rate  when  a  very  irregular  crack  front  is  present  or  where  there  is  elastic  and 
plastic  crack  tip  deformation,  which  in  these  particular  tests  is  significant. 

Since  the  effects  of  the  crack  face  closure  varied  among  specimens,  the  plan 
was  to  observe  the  potential  at  the  final  precrack  load,  and  then,  during  the  actual 
test,  assume  that  fracture  occurred  when  this  potential  was  exceeded.  Figure  A. 19 
shows  actual  signatures  obtained  from  two  tested  specimens  to  date,  namely  JC-31  and 
32.  The  basic  difference  between  the  tests  was  that  JC-31  was  tested  at  a  10-in. 
and  JC-32  at  a  12-in.  drop.  Note  that  their  general  appearance  is  similar;  there 
is  initially  a  very  large  negative  signal  (actually  overloading  the  recorder)  fol¬ 
lowed  by  a  positive  phase.  Based  on  the  previous  static  and  fatigue  work,  there  is 
no  reasonable  explanation  for  the  large  negative  peak.  In  Figure  A. 20  this  PD  trace 
is  shown,  along  with  the  corresponding  load  history  for  the  drop  tests  on  JC-32. 

The  approach  was  to  disregard  this  negative  peak  and,  based  on  fatigue  data  assume 
that  crack  initiation  essentially  occurs  when  the  trace  becomes  positive  (or  roughly 
5  mV).  This  occurs  at  about  5.8  ms,  at  about  the  time  that  the  applied  load  begins 
to  drop  off,  which  provides  some  indication  that  these  data  are  reasonable.  Note  in 
Figure  A. 19  that  for  specimen  JC-31  the  signal  becomes  positive  at  roughly  2.9  ms. 

The  applied  load  was  not  determined,  unfortunately,  due  to  faulty  circuitry.  However, 
specimen  JC-27  was  tested  at  similar  conditions  and  the  applied  load  for  that  test 
is  combined  with  the  PD  signature  in  Figure  A.21.  Again,  crack  initiation  is 
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PD  JC-31 


TIME  (ms) 


J _ L 
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Comparison  of  Potential  Difference  Record  with  Applied 
Load  for  Tests  at  Comparable  Conditions 


Indicated  when  the  signal  becomes  positive.  This  Indication  comes  early  In  the 
loading  phase  about  the  time  of  maximum  load.  It  could  not  be  certified  that  JC-31 
sustained  crack  extension  because  It  was  used  In  subsequent  tests  before  It  was 
broken  and  examined.  It  was  determined,  however,  that  specimen  JC-27  did  not  frac¬ 
ture.  In  the  static  and  fatigue  crack  work  this  procedure  looked  promising,  yet 
In  the  dynamic  tests  the  results  were  questionable.  There  Is  no  explanation  for 
the  large  signal  In  excess  of  1  V  obtained  during  the  tests  since  during  fatiguing 
the  signal  was  comparatively  smaller.  There  was  roughly  a  1/4-ln.  crack  extension 
and  the  voltage  change  was  only  0.23  V.  Also,  there  Is  no  reason  ^or  the  Initial 
large  negative  spike.  The  test  fixture  was  always  Insulated  from  ground,  yet  all 
the  Instrumentation  Is  common  to  ground  and  this  may  be  causing  some  problems.  The 
technique  needs  more  exploration  before  a  valid  assessment  can  be  made. 
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